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ABSTRACT 

We present new results from our study of the X-rayed outflow of the z = 3.91 gravitationally lensed 
broad absorption Hne (BAL) quasar APM 08279-1-5255. These results are based on spectral fits to all 
the long exposure observations of APM 08279-1-5255 using a new quasar-outflow model. This model 
is based on CL0UD\{^ simulations of a near-relativistic quasar outflow. 

The main conclusions from our multi-epoch spectral re-analysis of Chandra^ XMM-Newton and 
Suzaku observations of APM 08279-1-5255 are: 1) In every observation we confirm the presence of two 
strong features, one at rest-frame energies between 1-4 keV, and the other between 7-18 keV. 2) We 
confirm that the low-energy absorption (1-4 keV rest-frame) arises from a low-ionization absorber 
with log (A^n/cm"^) ~ 23 and the high-energy absorption (7—18 keV rest-frame) arises from highly 
ionized (3 < log ^ ^ 4; where ^ is the ionization parameter) iron in a near-relativistic outflowing 
wind. Assuming this interpretation, we find that the velocities on the outfiow could get up to ^ 0.7c. 
3) We confirm a correlation between the maximum outflow velocity and the photon index and find 
possible trends between the maximum outflow velocity and the X-ray luminosity, and between the 
total column density and the photon index. 

We performed calculations of the force multipliers of material illuminated by absorbed power laws 
and a Mathews-Ferland SED. We found that variations of the X-ray and UV parts of the SEDs and 
the presence of a moderate absorbing shield will produce important changes in the strength of the 
radiative driving force. These results support the observed trend found between the outflow velocity 
and X-ray photon index in APM 08279-1-5255. If this result is confirmed it will imply that radiation 
pressure is an important mechanism in producing quasar outfiows. 

Subject headings: cosmology: observations — X-rays: galaxies — galaxies: active — quasars: absorp- 
tion lines 



1. INTRODUCTION 

Th e existence of a Mbh ~ relation in nearby galaxies 
(e.g., Ferrarese & Merritt [2000) suggests that a feed- 
back mechanism exists regulating the co-evolution be- 
tween the massive black hole at the center of a galaxy 
and the formation of its bulge. Quasar outfiows could 
transport a fraction of the central black hole binding en- 
ergy to the host galaxy by possibly shocking against the 
host interstellar medium (ISM) and th erefore may be re- 
sponsible for the Mbh — relation (e.g., King|2011 ). Fast 
and powerful AGN jets and winds can be formed mainly 
by magnetic or radiative driving. Magnetically driven 
jets are thought to be responsible for the production of 
bubbles and cavities in galaxy clusters. The sources of 
these jets are found to be radio-loud (RL) AGNs located 
in the central giant elliptical galaxies (see 
Nulsen [2007 and references therein). Alt 



McNamara & 
rough there is 

compelli ng evidence for the exis tence of jets in RL AGNs 
(see e.g.. Bridle & Per ley 1984 and references therein), 
the presence of jets or magnetically-driven winds in ra- 
dio quiet AGNs is still debated. The mechanism driving 
winds in the majority of AGNs therefore is currently un- 
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''cloudy is a photoionization code designed to simulate conditions 
in interstellar matter under a broad range of co nditions. We hav e 
used version 08.00 of the code last described by Fcrland al.' (19981. 
Th e atomic database used by CLOUDY is described in l<'crg uson et] 
|al.| | |2001| and http: / /www. pa. uky.edu/ ~verner/atom.html, 



certain since the majority of AGN are radio quiet. ^ It is 
important to note, however, that a recent study claims 
observational evidence of magnetically-driven winds in 
high signal-to- noise ratio (S/N) spectr a of nearby AGNs 
(Fukumura et al. 2010[). Specifically, Fukumura et al. 
( ^10| ) find that a MHD model can explain the distri- 
bution of absorption as a function of ionization param- 
eter (dlog Nn/dS^) in the X-ray spectra of five nearby 
AGNs (mostly Seyfert 1 galaxies at z < 0.1). In is 
not clear if this result can be extended to high redshift 
sources or sources with relativistic winds where we expect 
strong and likely discontinuous gradients of the outfiow 
velocity and therefore we do not expect to find smooth 
dlog Nii/dS, distributions as required in their model. 

At moderate to high Eddington ratios (L/LEdd ?J 0.2), 
which may be common in high redshift sources z > 2 



(Kollmeier et al. 20061, we expect that radiative driv- 



ing become s more important than magnetic driving ( Ev- 



erett 



2005|) 



Recent observations of the high redshift 
traluminous infrared galaxy (ULIRG) SMM 
J1237-f6203 suggest a large-scal e powerful wind in thi s 
object driven by AGN radiation ( Alexander et al.|2010 ). 
These observations are suggestive that in high redshift 
AGNs (z > 2) radiatively driven winds may be more 
common than the jets found in radio galaxies. Radia- 
tive driving is also the favorite mechanism for explaining 
the formation of outflows observed in the UV from BAL 
quasars. Evidence of the later is found in the correla- 

3 RL AG Ns represent between 10-20% of the w hole AGN popu- 
lation (e.g., |Ivezi6, et al.|2002l [Jiang et al.|2007[ l 
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tion between changes in the maximum velocity of the 
outflow (C IV absorption lines) as a function of varia- 



tions in the so ft-to-hard spectral slope aox (e.g., Gib 



son et al. 



2009[ ) . The outflows observed in the UV from 

BAL quasars usually have moderate velocities (< O.fc), 
however, t here have been reported cases with velocities 



of ^ 0.2c ([Rodriguez-Hidalgo et al. |2007). Observations 
( [Hamann et al.|2008p of the quasar J f 05400.40-^034801.2 
imply wind velocities of ^ 0.8c based on measurements 
of the C IV line. We note, however, that even though this 
source reveals blue-shifted broadenings > 4000 km s~^ in 
C IV lines, it is not conside red a BAL quasar in th e strict 
sense of the definition (i.e. Weymann et al.|[l991 Sev- 
eral of the reported low S/JN cases of high v elocity AGN 
outflows have led to spurious detections (see Vaughan & 
Uttley [2008 and references therein). A handful of mod- 



erate S/IN X-ray spectra of BAL and mini-BAL quasars 
have been o btained that show strong signatures of fast 



winds (e.g., Chartas et al. 2002 2007a|b 2009a). The 



velocities of the outhowing X-ray absoroing material in 
some cases are found to be (> 0.7c) indicating a dif- 
ferent dynamical state than the observed outflow in the 
UV band. The fast variability and high ionization state 
of the winds observed in the X-rays could also be indi- 
cating that these X-ray outflows are originating closer to 
the central s ource than the outflows observed in the UV. 
Recently dSaez, Chartas & Brandt! |2009| [Chartas 



et al. 
tEe 



2009ap we presented observational evidence for 
presence of a powerful outflow in BAL quasar 



a powertui outnow m 
APM 08279-1-5255 based on the analysis of the X-ray 
spectra of this object. APM 08279+5255 is a high red- 
shift quasar which is unusually bright due to gravita- 
tional lensing, providing the exceptional opportunity to 
study it with high-quality spectral data. Our conclusions 
from these studies were that the outflow was accelerated 
to rclativistic speeds, and was releasing kinetic energy at 
a rate comparable to the bolometric luminosity of this 
source. The latter indicated that this outflow can re- 
lease an important fraction of the energy produced by 
black-hole accretion into the surrounding galaxy. In our 
study of APM 08279-1-5255 we also found evidence that 
the outflow could be radiatively accelerated by the cen- 
tral source. 

This paper is laid out as follows. In ^ we present 
the main conclusions extracted from recent stud- 
ies that focused on long X-ray observations (~100 
ks) of APM 08279+5255 performed with Chandra, 
XMM-Newton, and Suzaku. We also expand on the pre- 
vious studies in two directions. First in ^ we describe 
how the use of photoionization models may help us bet- 
ter constrain the physical properties of the outflowing 
X-ray absorbing gas. Second in ^ we discuss how the 
SED of the central source may influence the dynamics of 
the outflow. 

Unless stated otherwise, throughout this paper we 
use CGS units, the errors listed are at the 1-a level, 
and we adopt a flat A-dominated universe with Hq = 
70 km s-^Mpc-i, = 0.7, and VLm = 0.3. 



* The commonly used definition of BAL quasar (i.e. [Weymann] 
et al. |1991| does not consider outflow velocities > 2 5. OOP km s ^ . 
BA L wmds have also been found in the X-ray band fChartas et al. 
2002, 2007b), however, the X-ray faintness of BAL quasars (GaF 
lagher ct al., 2006(1 has limited the number of cases with observea 
A-ray BALs to a lew. 



Fig. 1. — Deconvolved image (left panel) of the 2002 February 
24 Chandra observation of APM 08279+5255. North is up, and 
east is to the left. Best fit PSF model (right panel) to the same 
observation of APM 08279+5255. 

2. THE FAST OUTFLOW OF APM 08279+5255 

This section contains information extracted from past 
studies of APM 08279+5255. In g]T] we concentrate 
on general information that is relevant to understanding 
the origin of the fast outflow of X-ray absorbing mate- 
rial observed in this object. In (S2.2), we describe the 



most important results extracted from eight deep X-ray 
observations of APM 08279+5255. 

2.1. Properties of the fast outflow of APM 08279+5255 

APM 08279+5255 is a gravitationally len sed BAL 
quasar at redshift z=3.911 (Downes et al. 1999). Its bolo- 



metric luminosity estimated as Lboi = 7 x 10^^ 
(hl is the lens amplification factor; Irwin et al 



1998 



Riechers et al. 2009) makes this source one of the most 



luminous in the universe^, despite the value of the lens 
amplification factor. The presence of high-ionization 
BAL features produced in C iv, O vi, N v and Ly a 
transitions indicates that it is a high ionization BAL 
quasar (HiBAL). The C iv BAL of APM 08279+5255 
contains multiple absorption troughs with widths rang- 
ing between 2000 km and 2500 km s^^. The maxi- 
mum observed outflow velocity of the C iv absorption is 
Umax(civ) ~ 0.04c relative to the source systemic redshift 
of z 



3.911 (Srianand & Petitjean 2000) 



Gravitational lensing of APM 08279+5255 produces 
three images with the two brightest ones separated by 
^ 0.5" (see Figure fl]). The time-delay between images 
A and B of APM 08279+5255 is estimated to be of the 



order of a few hours (e.g., Munoz, Kochanek & Keeton 
2001] ). The Eddington ratio of APM 08279+5255 can be 
estimated from an assumed value of F ~ 2 (where F is 
the X-ray photon index) , whic h is close to the mean value 
estim ated from observations (Saez, Chartas & Brandt 
|2009j iChartas et aL,2009a| , and using the /^boi/^Edd vs. 
I' correlation found in RQ qu asars (e.g., Wang et al.|2004 ' 

1' cor- 
0.2, 



|Sh emmer et al. 



relation tor AP 



2006 



. . 2008|). The Lboi/^Edd vs 

T5rD527g+5255 implies Lboi/^Edd 
and therefore A/bh ^ IO^^Ml^-^^o- A second indepen- 
dent method of estimating the Eddington ratio through 
the use of the observed width of the C iv line indi cates a 
black-hole mas s of A/bh ^ lO^^A^^^Af© (see e.g 



Riech- 



ers et al. 20091. The Eddington ratio derived using the 
second method is iboi/-^Edd ^ 2, implying the source is 
accreting at a super-Eddington rate. However caution 

^ As a reference t he most luminous qu asars in the SDSS have 
Lhol ~ 1O"L0 (e.g., [Ganguly et al.|2007 l. 
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TABLE 1 

Log of observations of APM 08279+5255. 



Date" OBS ID Telescope Instrument Exposure Net exp Not counts'' fo.s-s'^ 



2002-02-24 
2002-04-28 
2006-10-12 
2006-10-12 
2006-11-01 

2006- 11-01 

2007- 03-24 
2007-03-24 

2007-10-06 

2007- 10-22 

2008- 01-14 



(Epochl) 
{Epoch2) 
(OBSl) 
(OBSl) 
(OBS2) 
(OBS2) 
(OBS3) 
(OBS3) 
{Epoch3) 
(Epoch4) 
(Epoch5) 



2979 


Chandra 


ACIS BI 


91.9 ks 


88.8 ks 


5,627±75 


5.44±0.15 


0092800201 


XMM-Newton 


EPIC pn 


102.9 ks 


83.5 ks 


12,820±139 


5.51±0.15 


701057010 


Suzaku 


XIS FI 


102.3 ks 


71.3 ks 


6,071±90 


5.58±0.34 


701057010 


Suzaku 


XIS BI 


102.3 ks 


71.3 ks 


2,325±64 


5.39±0.43 


701057020 


Suzaku 


XIS FI 


102.3 ks 


67.9 ks 


5498±87 


5.18±0.32 


701057020 


Suzaku 


XIS BI 


102.3 ks 


67.9 ks 


2,181±62 


5.15±0.37 


701057030 


Suzaku 


XIS FI 


117.1 ks 


86.4 ks 


4750±80 


5.55±0.49 


701057030 


Suzaku 


XIS BI 


117.1 ks 


86.4 ks 


2,918±71 


5.63±0.39 


0502220201 


XMM-Newton 


EPIC pn 


89.6 ks 


56.4 ks 


11,400±114 


6.04±0.16 


0502220301 


XMM-Newton 


EPIC pn 


90.5 ks 


60.4 ks 


16,698±133 


8.06±0.15 


7684 


Chandra 


ACIS BI 


91.7 ks 


88.06 ks 


6,938±83 


5.89±0.21 



^ The date is also described by the term in parentheses. 

^ Background-subtracted source counts including events with energies within the 0.2—10 kcV band. See §2 of |Chartas et ah] 
j2009at and §2 of |Saez, Chartas &: Brandt] |2009| for details on source and background extraction regions used for measuring 

^ The fluxes (in units of 10""^"^ ergs cm~^ in the 0.5-8 kcV observcd-framc band arc obtained using the best-fit absorbed 

power law model with two absorption lines. The details of these fits arc presented in Table 3 of |Saez, Chartas fc Brandt| 
j2009[ and in Table 2 of [Chartas et ar!] l |2009a[ . The fluxes have been corrected for Galactic absorption. 



must be taken with this result because the estimates of 
Mbh based on C iv emission-hne widths is subject to 
a lot of systematic uncertainties (i.e., it contains non- 
virial components such as the d isk wind that can con - 
tribute to the broadening, e.g., IShemmer et al. 2008). 
There are also conflicting reported estimates of the mag- 
nification para meter (/xl)- Optic al observations indicate 
fiL 100 (e.g, Egami al. 2000), while observations o f 
CO emission yield /zl ~ 4 (e.g., Riechers et al. 2009). 
In this work we try, as far as possible, to present our 
results as a function of the magnification factor, however 
when t his is not possibl e a value of fiL = 100 is assumed 



(as in, Egami al. 2000). The absorption corrected op- 



tical to X-ray powe r-law slope^, of APM 08279-^5255 is 



1.8 (e.g., Dai et al. 2004). The relatively soft 



SED of APM 08279+5255 combined with its nearly large 
Eddington ratio provide ideal conditions f or the produc 



tion of radiatively driven winds (see e.g., Everett 2005 
and @. 

2.2. Results from eight deep X-ray observations of 
APM 08279-f-5255 

This section summarizes results from 2 Chan- 
dra, 3 XMM-Newton and 3 Suzaku observations of 
APM 08279+5255 (8 in total; see Table [T| for details). 
The data reduction of the Chandra and AMM-Newton 



observations is described in Chartas et al. ( 2009a I and 
the reduction of the Suzaku observations in |Saez, Char- 
tas & Brandt ( ,2009 ). Note that from here on when we 
refer to a particu lar observation we will use the same 
names defined in 



Saez, Chartas & Brandt (20091 



Chartas et al. (2009ap. I'herefore the three Suzak 



and 

ob- 

servations are referred to as OBSl, 0BS2 and 0BS3, 
the three XMM-Newton observations are referred to as 

® Oox is defined as the slope of a hypothetical power law ex- 
tending between 2500A and 2 keV in the AGN rest frame, i.e. 



: log 



F^(2keV) 



. 0.3838 log 1^(2^ 



u{2keV) 

■F„(2500A) / '°S^{2500A) — "-^ "-J" f „ f250nA1 

The Oox— ^uv relation described in jStratevaet al.|2005[|Steffen| 
|et al ||2006|l predicts t hat Onx = ( -0 137 ± 0.008) log i^F^nnA + 
\:i.b38 ± U.240). Using [irwin et al.] | [T99 8> work we estimate that 
log Z2500A ~ 3l.7 (/ii ~ lUO) and therefo re we obtain Oqx = —1.7+ 
0.3 in accordance with |Dai et al.| | |2004[ | value. 



Epoch 2, Epoch 3 and Epoch 4 and the two Chandra 
observations are referred to as Epoch 1 and Epoch 5 (see 
Table [l] for details). 

Our analysis of the X-ray observations of the BAL 
quasar APM 08279+5255 indicates strong and broad ab- 
sorption at rest-frame energies of 1-4 keV (low-energy) 
and 7-18 keV (high-energy). The medium producing 
the low-energy absorption is a nearly neutral absorber 
with a column density of log (iVn/cm"^) ~ 23. However, 
since APM 08279+5255 is at a high redshift of z = 3.91 
it is difficult to constrain the ionization state of the low- 
energy absorber. The absorption signatures of this com- 
plex low energy absorber are shifted outside the range of 
the observed X-ray band (typically between 0.3-10 keV). 
The high-energy absorption features are easily detected 
in the 8 deep X-ray observations of APM 08279+5255. 
In every case these features sa, t isfy EW/gF,w ^. 3 (se e , 



Saez, Chartas & Brandt 



2009 



Chartas et al. 



2009a) 



and theref ore they satisfy realistic li mits tor signihcant 
detections ( |Vaughan fc Uttley |[2008l ) . 

Chartas el al. 2002 have interpreted the high-energy 
X-ray BALs as being produced by absorption of highly 
ionized iron such as Fe xxv Ka (Is^ — ls2p; 6.70 keV) 
and/or Fe xxvi (Is — 2p; 6.97 keV) launched very near 
an ionizing compact central source. Evidence to sup- 
port this interpretation is the significant variability in the 
strength and energy of the X-ray BALs over short time- 
scales present in many of the observations. This variabil- 
ity is found to be as short as ^ 3 days rest-frame in the 
Chandra and XMM-Newton observations, and ~ month 
in the Suzaku observations. This fast variability presents 
a strong argument in favor of a wind that originates from 
a distance of a few times the Schwarzschild radius^. 

* EW is the equivalent width obtained using the best fitted 
model of the high energy absorption feature. The equivalent width 
(EW) is defined as EW = f "p ^ dE, where Fc is the continuum 
flux and Fe is the flux in the absorber. 

^ The time-scale of the flux variability provides an indication 
of the outflow launching radiu s from the central engine (see e.g., 
baez, Chartas k. Brandt||2009| . The size of this region should be 
approximately the time scale of the variability times the speed of 
light. Variability over a period of a ^week implies a launching ra- 
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TABLE 2 

The minimum and maximum energies and velocities of the high-energy absorption 
FEATURES IN APM 08279+5255". 



OBS 


Instrument 


-^min 




■^min 


^max 


r 






[kcV] 


[keV] 


[c] 






OBSl 


XIS FI 


7.52±0.36 


12.16±0.52 


0.11±0.05 


0.53±0.03 


1.94±0.04 


OBS2 


XIS FI 


7.88±0.48 


11.57±0.51 


0.16±0.06 


0.50±0.03 


2.02±0.03 


OBS3 


XIS FI 


7.51±0.32 


11.70±1.06 


0.11±0.04 


0.51±0.07 


1.94±0.05 


OBSl 


XIS BI 


7.24±0.51 


12.47±1.12 


0.08±0.07 


0.55±0.06 


1.95±0.07 


OBS2 


XIS BI 


8.74±0.58 


12.24±1.16 


0.26±0.06 


0.54±0.07 


1.91±0.06 


OBS3 


XIS BI 


7.53±0.82 


12.40±0.92 


0.12±0.11 


0.54±0.05 


1.92±0.05 


Epochl 


ACIS S3 


8.05±0.11 


10.64±0.19 


0.18±0.01 


0.43±0.01 


1.74±0.03 


Epoch2 


EPIC pn 


7.54±0.42 


15.04±0.97 


0.12±0.06 


0.67±0.04 


1.89±0.03 


Epochs 


EPIC pn 


6.42±0.48 


17.94±1.81 


< 0.03 


0.76±0.04 


2.03±0.04 


Epoch4 


EPIC pn 


6.94±0.20 


16.31±0.90 


0.04±0.03 


0.71±0.03 


2.11±0.02 


Epoch5 


ACIS S3 


7.36±0.30 


15.35±1.13 


0.09±0.04 


0.68±0.04 


1.94±0.04 



The model used to estimate Emin and Smax is a power-law with Galactic absorp- 
tion, intrinsic absorption, and t wo absorption lines (APL+2AL ). The fits of this 
model are presented in Table 3 of Saez, Chartas & Brandt ( 2009 1 and in Table 2 of 
Chartas et all (|2009a|). 



In the Suzaku observations the variabiUty of the 
high energy absorber is present at energies close to 
7 keV rest-frame. On the other hand the Chandra and 
XMM-Newton observations show variabihty of the high- 
energy absorption profiles at both: low energies (i.e. 
rest-frame energies ~ 7 keV) and high energies (i.e. 
rest-frame energies > 9 keV). A time- variable outflow 
provides a plausible explanation for the changes in shape 
on the absorption fea tures in past X-ray observations of 
APM 08279-H5255 (Chartas et al.||2002[ |Hasinger et~ar 
2002 1. An even stronger case is presented in Chartas et 
arn [2009a[ ) and [Saez, Chartas fc Brand"t] ( |2009p , where 
it can be seen clearly that in some cases the appear- 
ance of the high-energy absorption feature can take the 
shape of a notch, an edge or two absorption lines. Our 
spectral analysis of the 8 deep observations also indicates 
variability of the parameters defining the high energy at- 
tenuation which is likely due to a change in the outflow 



Chartas et al. 



velocity of the absorb er (Saez, Chartas & Brandt 2009 
Ihe short time-scale ( 



2009a[ ) . T'he short time-scale ('^week m 
the rest-frame) of the variability combined with the high 
ionization of the absorbing material which is moving at 
relativistic speeds imply that the absorbers are launched 
from distances < 10 R s from the central source (see e.g 
§4 of Saez, Chartas & Brandt 200 9|. The short time-scale 
of this variability in the S uzaku (jSaez, Chartas fc: Brandt 
'2009| and XMM-Newton ( [Chartaset al.|2009ap observa- 
tions also indicates that this absorber should be strongly 
accelerated. The global covering factor of these winds 
should be low < 20% based on the absence of emission 



features from the outflowing ionized gas (see e.g. Chartas 
et al.||2009a|). 
TFe 



minimum 



and maximum projected velocities 
) of the outflow are estimated from the mini- 



mini ^max) the 



mum and maximum energy ranges (£", 
high-energy absorption features in APM 08279-1-5255. 
We obtained -Emin and E'max from our spectral flts as- 
suming first the two absorption-line (APL-I-2AL) model. 
Specifically, based on the best-fit values of an absorbed 

dius ^ 6-Rs. In the last expression we estimated the Schwarzschild 
radius (Rg) of APM 08279-1-5255 using A/bh ~ W^°Mq (obtained 
assuming /^^ ~ 100; see |2.1[ l. Possible flux variability caused by 
the time-delay between the images is shorter than -Rg/c. 



power-law model with two absorption lines (from Table 3 
of Saez, Chartas & Brandt 2009 and Table 2 of Char- 
tas et al. 2009), we obtain £^mi n = -£^absi ~ 2gabsi an d 

As 



abs2 



2a. 



abs2- 



Chartas et al. (2009a) 



we estimate the line of sight projected velocities" timin 
and Wmax assuming the absorption arises from highly 
blueshifted Fe xxv Ka (Ei^b = 6.7 keV). In Tabled we 
show the values obtained for -Emim ^-max? ''^min and Umax- 
In this table we also add a column with the fitted values 
of F (based on the APL-I-2AL model). Using -Emax ~ 
18 keV (the maximum observed value), we constrain the 
maximum angle^*^ between our line of sight and the wind 

direction to be less th an 22°. 

As we indicated in Chartas et al. ( 2009a[ ) (and from 
Table [2]) there is a hint that changes in the photon in- 
dex (Ij may be positively correlated with the changes of 
the maximum velocity of the outflow (i^max)- This possi- 
ble trend between F versus Vmax is shown in Figure 10 of 

3]we recalculate the velocities 
el based on cloudy simula- 



Chartas et al 
found in 



TSBTe 



(I2009a[). In 
iMwIth a mo 

tions of a near-relativistic outflow. In Q we provide an 
interpretation of the possible trend between maximum 
outflow velocity and F. 

3. PHOTOIONIZATION MODELS OF NEAR-RELATIVISTIC 
OUTFLOWS 



3.1. Motivation 

From the X-ray ana lysis of our 
APM 08279-H5255 (see 



8 observations of 



Saez, Chartas & Brandt 2009 



Chartas et al. 2009a for more details), and pa st obser 



vat ions of th is object and PC 1115-1-080 (e.g., Chartas 
et al. 2007a), it became clear that a more sophisticated 
spectral model was needed to fit the X-ray BAFs. Cur- 
rently there are no proper tools available in the spectral 
fitting package XSPEC to model absorpti on profiles result- 
ing from near -relativistic outfiows. In Saez, C hartas fc| 
Brandt (2009) we used photoionization models created 
from the photoionization code xstar (see, e.g., |Kall-| 

The Doppler-shift formula predicts that given a fixed ratio 
of -Eiab/-Bobs=filo (where -Eiab and -Bobs is the energy of the ab- 
sorption line in the rest-frame and observed-frame respectively) the 
maximum angle between our line of sight and the wind direction 

is given by Smax = 008^1(^/1^^/?^). 
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man & 



Bautista|2001[ ) to fit tfie broad absorption feature 
)etween7-18keV rest-frame in APM 08279. This 



found between ' 
spectral analysis did not include relativistic corrections, 
and does not contain a realistic model of the outflow. In 
this section we describe new software code that we have 
developed that provides a more realistic description of 
radially accelerated near-relativistic outflows. Our main 
goal is to better constrain important parameters describ- 
ing near-relativistic outflows; among these are, the ve- 
locity profile, the ionization parameter and the column 
density of the observed X-ray wind. 

3.2. Description of the code 

Our quasar outflow code is based on a multilayer 
approach which mimics the absorption and scattering 
through a near-relativi stic outflow using the photoion- 

An existing 



ization code CLOUDY (Ferland al. 



1998) 



quasar outflow code called XSCORT developed by Schurch 
& Done 2007 follows a similar approach and is based on 
the photoionization code XSTAR. 

In our quasar outflow code we use cloudy simulations 
to approximate the absorption signature produced by an 
outflowing medium that is radiatively accelerated from a 
central source. We describe the outflowing medium with 
a set of absorption layers with specific velocity profile 
and ionization states. The details of the approximations 
and assumptions used in our quasar outfiow model are 
described in the following paragraphs. 

3.2.1. The case of one absorption layer 

We approximate the attenuation through the outflow- 
ing ionized gas as the absorption through a series of layers 
of different velocities and ionization states. We begin by 
calculating the output absorption profile assuming the 
wind is made up of a single layer and then generalize 
our model by dividing the outflowing wind into multiple 
layers. We assume the layer is at a distance R from the 
continuum source, has a thickness AR <C i?, a density n, 
and is moving at a velocity v away from it. We further 
assume that spectral energy distribution of the source is 
a power-law (L^ oc u") with spectral index a = —1.0. 
From here on, unless mentioned otherwise, this is the 
SED that we wi ll assume in our simulations. Th e ioniza- 
tion parameter ( Tarter, Tucker fc Salpeter|1969 see Ap- 
pendix IA| measured in the absorption layer's rest-frame 



IS 



(1) 



where primed quantities n' in' = ^n) and Fj refer to 
the density and the incident ionizing flux in the layer's 
rest-frame, respectively. Non primed parameters are as- 
sumed to be in the rest-frame of the luminous source. 
Using the fact that lyjv^ is a Lorentz invariant, the in- 
cident ionization flux in a layer's rest-frame is given by^^ 

Since oc I^, u/u' = D; where D = 'y^^ (l—v / c)~^ and Ii, = 
7,yi/3/jy'3 = /^,Z)3. Therefore F'j = /J/ F^.dv' = D'* f^^ F^du, 

where i/g = 1 Ry and v'^ = 1000 Ry; vq^\ = Du'q j. Conversely, 
to calculate the incident flux in the layer's rest-frame, we integrate 
the flux in the luminous source rest-frame, and multiply by a factor 
of D~^. Notice that if the flux is a power law, i.e. F,^ = av°' , 



F'j^^\l~ v/cfFi 



(2) 



The column density (in any reference frame^^) of the 
layer is given by A'h — nAr. The incident flux on a 
layer is calculated between 1-1000 Ry in the rest-frame 
of the moving layer. ^'^ From the incident flux in the rest- 
frame of the layer and from the use of equation ([T]) we 
calculate the ionization state of this layer. We next mput 
the derived ionization state and column density of the 
layer into the photoionization code CLOUDY to obtain 
the output spectrum of the layer in its rest-frame and 
transform it to the rest-frame of the luminous source. 

3.2.2. A model to describe X-ray absorption profiles. 

The X-ray spectra that we have obtained from obser- 
vations of several BAL quasars do not sho w any emis- 
sio n lines associated with the X-ray BALs (Chartas et 



al. 2009a I. This is expected since BAL quasar outflows 



then F'j = aD~Hn{u[/iy[,) if a 



-1, and Ft 



m general have relati vely small global covering factors 
( Hewett fc Foltz|2003 )■ Therefore, we will concentrate on 
absorption-dominated outflows neglecting emission fea- 
tures in our model. For our cloudy code, we will as- 
sume an approximately constant ionization parameter in 
the outfiow. 

We assume that the region of the outflow that con- 
tributes to the observed Fe xxv and/or Fe xxvi absorp- 
tion is localized in a relatively thin multilayer of thickness 
AR at a distance R from the black hole. The justifica- 
tion behind the assumption that AR <C i? in our model 
is that the presence of Fe XXV is likely to occur in a rel- 
atively narrow region away from the source where the 
ionization parameter and density are such that Fe xxv 
is abundant at these distances. Even though the ioniza- 
tion parameter is expected to fall off with distance along 
the fiow as 1 /i?^ it remains relatively constant within the 
multilayer since AR <C i?.^'^ To compensate for changes 
in flux across the multilayer due to relativistic effects 
(e.g., beaming effects) we adjust the density Ui of each 
layer such that the ionization parameter remains approx- 
imately constant within the multilayer AR. In the case 
that the region of the outflow that we are modeling is 
not physically thin, our approach will still be valid as 
long as the ionization parameter along this region does 
not change appreciably. 

From equation (n]) the ionization parameter of each 
layer is = AttFIJti^; where Fj ^ is the incident flux in 
the rest-frame of layer i estimated from the flux coming 

!/o"+^)/{a + 1) if o ^ -1. Therefore in the case that the SED is 
a power-law, the ionizing flux in the rest-frame of the absorbing 
layer can be found by integrating the flux in the luminous source 
rest-frame between 1 Ry and 1000 Ry, and then multiply by D"~^ 
where a is the spectral index. 

12 iVn = nAr = (n7)(Ar/7) = n'Ar' 

1^ The Rydberg is a unit of energy deflned in terms of the 
ground-state energy of an electron for the hydrogen atom, 1 Ry 
13.6eV. 

1* In our simulations we assume that the ionization luminos- 
ity is Lj = 10** erg s~^, therefore for a layer with log ^ = 3.0 . 
n = 10^2 cm"'' (initial density assumed in our runs, see ^3.2.4 1, 
and log (Afjj/cm^^) = 23, we will typically have that the ratio 
of thickness to distance of the layer is AR/R Ri 0.001 and conse- 
quently A^ —2. This means that the ionization parameter with 
typical value of log f ~ 3 varies through the outflow ~ 0.001 dcx 
due to the radial dependency. 
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from layer i — 1 after taking into account relativistic ef- 
fects. The density of each layer, is determined from an 
a priori defined velocity profile v{r). In our multilayer 
approach, the density of each layer is Ui cx 7,f (1 — Vi/c)^ 
for i = 0, l..{Ni — 1), were Ni is the number of layers. We 
note that the assumed density profile adjusts the ioniza- 
tion parameter so that it does not vary across the mul- 
tilayer due to relativistic effects. However, the ioniza- 
tion parameter will still decrease betwe en 0.02 — 0.40 dex 
mainly due to attenuation (see 5 3.2.3). For our case we 
concentrate on simulated outflows with log ^ > 3. As 
long as each layer^^ used in our calculations has a col- 
umn density of log (AA^Hi/cm~^) < 23 it can be assumed 
to be almost transparent (or th in) to radiation (see Ap- 
pendix |b|). As we describe in 5 3.2.3 the assumption of 



optically thin layers will allow us to more greatly reduce 
the time needed to compute the absorption signature of 
the multilayered outflow. 

In this work we assume that the outflow velocity profile 
will have a p-type form given by 



+ {v„ 



V^ninKNii/Nnf, (3) 



where A^h is the total column density of the simulated 
wind (0 < A'h < A^h)- We note that given the assump- 
tions described here, and the moderate S/N of the ab- 
sorption profile in X-ray, we do not expect to obtain 
significant constraints on the acceleration mechanism of 
the outflow by using this profile. 

3.2.3. Passing a continuum spectrum from one absorbing 
layer to the next 

When passing the continuum spectrum from one layer 
(i-1) to the next (i) we first transform the source 
rest-frame spectrum coming from layer i-1 to the rest- 
frame of layer i. In the rest-frame of layer i we remove 
absorption and add emission (in cases where we do con- 
sider emission) to the incident spectrum given the col- 
umn density and ionization parameter of the layer. The 
column density of layer i is obtained from the velocity 
profile, and its ionization parameter is calculated from 
the input spectrum in its rest-frame. Finally we trans- 
form the output spectrum coming from layer i back to 
the source rest-frame. 

Since the density of each layer is chosen to correct 
the attenuation of the flux due to relativistic effects, 
the ionization parameter will mostly decrease due to 
the absorption of the outflow across the layer. The de- 
crease of the ionization through the multilayered out- 
flow due to absorption will be greater in winds that are 
less ionized and with higher column densities. For ex- 
ample, a wind with log (A^n/cm"^) — 22.75 with ini- 
tial ionization parameter of log — 3.25 will have a 
final ionization parameter of log ^ 3.21, however for 
log (A^n/cm-^) = 22.75 and log = 3.75 the final ion- 
ization parameter is log f j ~ 3.73. Additionally a wind 
with log (A^H/cm-2) 23.75 and log = 3.25 will have 
log - 2.84, however if log (A^n/cm-^) = 23.75 and 
log = 3.75 then log ~ 3.53. 

We use 100 layers to generate the absorption profile of 
a simulated wind, therefore each layer has a column density of 
log ANiii < 22. 



After correcting the spectrum incident on layer i for 
the effects of relativistic beaming, we estimate the ab- 
sorption and emission from layer i using simulations per- 
formed with the photoionization code cloudy. These 
CLOUDY simulations are implemented using the unab- 
sorbed source spectrum (power-law with a = — 1 in our 
case) and various values of the ionization parameter. 
The main advantage of this approximation is a dramatic 
shortening in the time it takes to obtain a profile. This 
speed-up is possible because the photoionization runs at 
each layer use a library of preexisting cloudy runs (see 
next paragraph for details). If the flux emitted by the 
source is a power-law, the flux received by a moving layer 
is also a power-law with the same spectral index as the 
source as long as the SED is not significantly attenuated 
in an energy dependent way. In addition, if each layer 
of the multilayered media is approximately optically thin 
then the state of the gas depends primarily on the ion 



ization parameter and secondarily on the SED (Tarter, 
Tucker &: Salpeter|1969) . This dependence indicates that 



our approach is at hrst order a good approximation as 
long as the SED along the multilayered outflow does 
not differ appreciably from the incident SED at the first 
layer. However, for low ionization states (log ^ < 3) and 
large total column densities (log (A^n/cm"^) > 23.5), the 
SED along the multilayered outflow will deviate appre- 
ciably from the initial incident SED, and therefore, our 
approximation breaks down. 

In order to avoid invoking cloudy too many times 
(~ 100) to simulate the signature of the multilayered ab- 
sorber, we created a library of spectra to speed-up our 
calculations. This library contains a set of simulated 
transmitted spectra through absorbers having a range of 
ionization parameters between 1.8 — 4.1 dex (log ^ steps 
of 0.01 dex), column density of log (A^n/cm"^) = 19, 
density of nn = 10^^ cm~'^ and an input SED that has 
the form of a power law {L^ oc v") in the energy range 
between 1 — 10"* Rys with spectral index a = —1. The 
value of the density of the absorber (nn ~ 10^^ cm^'^) is 
chosen to be close to the value of the density predicted by 
^""-^ Dunicrical s imulations of radiative winds by Proga et| 
al. (2000 2004) ) . The simulations used to construct this 
TTBrary (from here on referred to as library runs), as well 
as all the simulated profiles described in the next sec- 
tions, are calculated using standard solar metallicities. 
We primarily used cloudy to perform the photoioniza- 
tion calculations but we also compared our results with 
those generated by the photoionization code xstar. 

In order to reduce the number of realizations in our 
library, we performed these runs assuming a specific tur- 
bulent velocity of Wturb = 1500 km s~^. Our choice of a 
specific turbulent velocity requires that the relative speed 
Av[ between layers i and j -I- 1 to be the same for every 
layer in the simulations of multilayered outflows^^ (i.e. 
Aw,- = Av'; see Appendix [C|). In our simulations we use a 
finite number of layers to describe the outflow, each layer 
having a slightly different velocity following the adopted 
velocity profile. The distribution of the velocity differ- 
ences between adjacent layers is assumed to be uniform 
and to extend between and a maximum velocity differ- 
ence of Aw'. The velocity difference of Av' is kept lower 

^® The relative velocities of two adjacent layers according to 
special relativity is Av'^ = (fi+i — Vi)/(1 — ViVi^i/c^). 
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than Wturb\/l2; where Wturb = 1500 km is th e as- 
sumed turbulent velocity (see Schurch fc Done|2007 ) . By 
selecting a turbulent velocity we set a lower limit on the 
number of layers used in simulating a BAL spectrum. 
The simulations making up our library are all performed 
using the same column density (log (iYn/cm^^) — 19). 
Consequently, the opacity of a layer i (i.e., the layer's 
rest-frame absorption) is calculated by first obtaining 
the opacity from the library at the ionization state of 
the layer, and second, by scaling the opacity obtained 
(from the library) to the column density of the layer^®. 
This approximation is justified given that in the range of 
densities of our simulations (8 < log nn ^ 12) the depen- 
dency of the opacity on is negligible, and therefore, 
the opacity of a layer i can be found from the library 
just by knowing its ionization state and column density. 
Emission spectra may also be obtained from our library 
runs, however for our study we mostly concentrated on 
the absorption spectra. We focus on the absorption spec- 
tra since the main goal of our study is to fit the spectral 
signatures produced by winds in which emission lines are 
relatively weak implying low global covering factors. Ad- 
ditionally, we assume that the line-of-sight covering fac- 
tor is equal to one. 
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Fig. 2. — X-ray spectra derived from multilayered CLOUDY simu- 
lations of near-relativistic outflows. The assumed input spectrum 
is a power-law with F = 2.0. The outflow has been accelerated 
between 0.1c — 0.7c using a p-type velocity proflle with p = 1.0. 
The total column density of the outflow is log (Afjj/cm"^) = 23.75. 
We have calculated the output spectra for 3 different values of the 
ionization parameter, log ^ = 3.0, 3.5 and 4.0. 



3.2.4. Testing our relativistic outflow code 

We first generate the absorption profiles for a linearly 
accelerated outfiow with Vmm ~ 0, Wmax = 0.3c [v oc .^^h), 
= 3 X 10^^ cm^^, initial density tiq = 10"^^ cm~^, 
and values of log ^ = 2.75, 3.00, 3.25, 3.50, 3.75 and 
4.00. These profiles have been generated to c ompare the 
result s of ou r code with those presented in [Schurch "SE] 
Done (2007). In order to make a fai r comparison with 
the results of Schurch & Done (2007), we use a library 
of XSTAR runs using a power-law model with a — —1.4. 
This library of XSTAR runs, has the sa me ch aracteristics 
as the CLOUDY library described in J 3.2.3 Our simu- 
lated X-ray absorption profiles are very similar to the 
ones computed by Schurch & Done 2007 under the as- 
sumption of a multilayer outfiow with spherical symme- 
try (see Appendix [d|). 

We generate a second run of simulations based on our 
results presented in §2.2[ For this purpose, we simu- 
late absorption profiles with i;„iin = 0.1c, Wmax = 0.7c, 
no = 1012 cm-3^ log (TYn/cm-^) = 23.75, and three 
values of the ionization parameter log ^ = 3.0,3.5,4.0. 
The outflow parameters have been chosen approximately 
in accordance to the estimated values based on one 



of the XMM-Newton observations of APM 08279 
(Epoch 3 



from 



2009a I 



5255 
To generate 
= 1 for these 



Chartas et al 

notch- type absorption profiles we chose p 
simulations. The resulting profiles are presented in Fig- 
ure [2] In Figure [3] we zoom in on the iron blend fea- 
ture which shows up at energies between 7.5—16.5 keV. 
The resonance absorption lines responsible for the fea- 
tures are mainly Fe xxv Ka (Is^ — ls2p; 6.70 keV) and 

For example if the velocity profile has i>min = 
and Umax = 0.7c we require that the number of layers 

Ni > (-Umax - fmin)/(^2-!;turb), 1-6. with Dtu^b = 1500 kms'^, 

the minimum number of layers that we can use for this case is 
forty. 

For thin layers, the opacity of a layer (t2) can be obtained 
by scaling the opacity of a different layer (ri ) by the ratio of the 
column densities of the two layers (i.e. T2 = TiAfH2/^Hi)- 
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Fig. 3. — X-ray spectra derived from multilayered CLOUDY simu- 
lations of near-relativistic outflows. This is a zooming of Figure[2]in 
the region where we find iron absorption. The assumed input spec- 
trum is a power-law with F = 2.0. The outfiow has been accelerated 
between 0.1c— 0.7c using a p-type velocity profile with p=1.0. The 
total column density of the outflow is log (Afjj/cm" ) = 23.75. 
We have calculated the output spectra for 3 different values of the 
ionization parameter, log ^ = 3.0, 3.5 and 4.0. For reference the 
observed blueshifted energies of the Fe xxv Ka and Fe xxvi Ka 
lines with v = 0.1c and of the Fe xxv Ka and Fe xxvi Ka Unes 
with V = 0.7c are highlighted. 

Fe XXVI Ka (Is — 2p; 6.97 keV). From Figures [2] and 
[3] it is clear that the strength of the iron blena fea- 
ture is larger for log ^ = 3.5. The strength of this fea- 
ture is large for ionization parameters in the range of 
3 ^ log C ^ 4, since at lower ionization levels the abun- 
dance of highly ionized iron atoms (e.g., Fe XXV and 
Fe xxvi) is too small to produce any appreciable fea- 
ture in the spectra. At log ^ ^ 4 the atoms in the gas 
become stripped of almost all their electrons, and there- 
fore, the resulting spectra contain very few bound-bound 
and bound-free absorption features. When log C ^ 3 
there is substantial absorption present at rest-frame en- 
ergies < 2 keV. As a consequence, in a heavily absorbed 
source like APM 08279-f 5255, the low-energy (< 2 keV) 
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Fig. 4. — Outflow velocity as a function of column density 
through the outflow for 3 different values of the the exponent of 
the p-type velocity proflle (see equation [sjl , p = 0.5 (dotted line), 
1.0 (solid line) and 2.0 (dashed line). We assumed the outflow has 
a launch velocity of 0.1 c, a terminal velocity of 0.7 c and a total 
column density of log (A''jj/cm~^) = 23.75. 
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Fig. 5. — Upper panel: X-ray spectra produced from multilay- 
ered CLOUDY simulations of near-relativistic outflows. The assumed 
input spectrum is a power-law with F = 2.0 (solid line). The out- 
flow has been accelerated between 0.1c— 0.7c using an exponent of 
the velocity profile with p=1.0. The total column density of the 
outflow is log (A'jj/cm"^) = 23.75. We have calculated the output 
spectra for 3 different values of the exponent p, p = 0.5 (dotted 
line), 1.0 (solid line) and 2.0 (dashed line). Lower panel; X-ray 
spectra resulting from the best fitted models used in Chartas et 
[al| (2009aj to fit the high energy absorption of APM 08279+5255 
m lipoch 3. The solid line is for the notch model, the dotted line 
corresponds to the two absorption line model and the dashed line 
is for the absorption edge model. 

absorbed component of the spectrum, possibly due to 
non-outflowing material, will overlap w ith the absorbed 
component produced by the wind ( e.g., Saez, Chartas &] 
Brandt||200"9| [Chartas et al.||2009a[ ). However, the high- 
energy absorption iron trough (above ~ 7 keV) is not 
contaminated by absorption lines from non-outflowing 
material, and consequently, it provides a clean signature 
of the outflowing material. 

The absorption profiles described in §2. 2| contain a di- 
versity of shapes, sometimes resembling notches, edges or 
two absorption lines. Given the variety of shapes of the 



absorption profiles we attempted to emulate these dif- 
ferent shapes by generating simulations of outflows with 
Vmin = O.lc and Wmax = 0.7c, log (iVn/cm-^) = 23.75, 
log ^ = 3.5 and three values oi p, p = 0.5, 1 and 2. The 
output spectra of the resulting velocity profiles (see Fig- 
ure |4] ) are presented in Figure [5j In the lower panel of 
Figure [5] we also show as a reference three model spectra 
used to fit Epoch 3. The three different models used are 
a notch, an edge and a two absor ption line model (see 
Table 2 fromlChartas et al.||2009a|). 



3.3. Fits to the observed spectra of APM 08279+ 5255 
using our quasar outflow model 

In this section we present results from fits to spectra 
obtained from XMM-Newton, Chandra and Suzaku ob- 
servations of APM 08279-1-5255 using our thin multilay- 
ered outfiow model. We created xspec table models^^ 
from a large number of simulated quasar absorption spec- 
tra produced from running our multilayered model over 
a large range of parameter space of the outflow param- 
eters. Specifically, the parameter space of our simula- 
tions covered five important properties of the multilay- 
ered outfiow model: minimum velocity of the outflowing 
gas Wniini maximum velocity Umax: ionization parameter 
^, total column density A^h and the value p which de- 
fines the outflow velocity proflle (see equation 3]). For 
our simulations we assume (unless mentioned otherwise) 
an initial density of the outflowing gas in the multilayer 
of no = 10^^ cm~'^. We note that our results are not 
sensitive to the assumed initial density of the outflowing 
gas. A change of the density by up to two orders of mag- 
nitudes from our assumed value does not result in any 
significant variation of our results. 

3.3.1. XSPEC table models for quasar outflows 

In order to fit the variety of absorption spectra of 
APM 08279-1-5255 we generated three different types of 
XSPEC table models that mainly differ in the parameter 
space of the outfiow properties that are covered. 

The first table model (windfull.tab) assumes the 
presence of only one outflowing component. We define 
as component i of an outflow a multilayer absorber at 
a distance R{ from the black hole, with a thickness Ai?i 
{ARi/R[ ^ 1), a velocity gradient of Avi across it, a col- 
umn density of A'h, and an ionization parameter of log ^j. 
We produced the table model windfull.tab by simulat- 
ing transmitted spectra through a range of outflows with 
minimum velocities (wmin) in the range of 0.0— 0.3c (a 
simulation for every 0.02c interval), with maximum ve- 
locities (vmax) in the range of 0.3— 0.9c (step intervals of 
0.04c), log (A^n/cm"^) in the range 22—24 (step intervals 
of 0.25 dex), log ^ in the range 2.8—4.0 (step intervals of 
0.2 dex), and the index p having values of 0.5, 1.0, 2.0 
and 5.0. 

In several observations, the X-ray BALs produced by 
highly ionized iron lines (above 7 keV rest-frame) appear 
to consist of two broad absorption components. The 
component covering low energies (< 9 keV rest-frame) 

A table model consist of a file that contains an A^-dimensional 
grid of model spectra with each point on the grid having been cal- 
culated for particular values of the A'' parameters in the model. 
XSPEC will interpolate on the grid to get the spectrum for the pa- 
rameter values required at that point in the fit. 
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TABLE 3 

Results from spectral fits to the three Suzaku observations of APM 08279+5255. 



Model" 



Parameter 



FI SPECTRUM* 
Values OBSl Values OBS2 Values OBS3 



BI SPECTRUM'' 
Values OBSl Values OBS2 Values OBS3 



MODELl. 



MODEL2., 



log (AfHi/cm-2) 
log {AfH2/cm-2) 
log ?2 

^max / C 

p 

log L2-10 
V 

log 

log (AfH2i/cm-2) 
log (AfH22/cm"2) 

log ?2 
^(1) /c 

min' 

j;''' /c 

(2) / 
Ic 

•^max/ 

log L2-10 



1 96+°-'''' 
22 gi+0-09 

23.54l«-i; 
o 44+0.16 

14+002 

'J-0'-'-0.02 
1 6+0-4 

46.92+0.08 
60.5/67 
0.70 



j+0.05 
-0.05 
,+0.16 
-0.20 
,+0.15 
-0.15 
-,+0.20 
-0.18 
j+0.28 
-0.18 
,+0.02 
-0.02 
j+0.03 
-0.03 
3+0.06 
-0.06 
3+0.03 
-0.02 

46.95+0.08 



1.98: 
22.82:1 
23.23j 
23.26j 
3.28;! 
0.12^ 
0.28:1 

0.39:! 

0.49"* 



55.1/65 
0.81 



+0.06 
-0.05 

22 q5+0--'-2 

00 00+0. 18 
ZO.OZ_(j 

o 20+0.10 
•""■^^-O.IO 
3+0.02 
0.02 

+0.04 



2.03 



0.19 



1.95 



D ^I ^ 
'^■'^-^-0.04 

46.91+0.07 



68.8/60 
0.21 



3+0.04 
-0.04 
7+0.07 
-0.07 

23.l0i0;?« 
22.98l°-|6 
3+0.12 

-0.16 
-1+0.03 
-0.03 

29+0 04 
'-'■^^-0.04 
41+0-05 
'-'-^-^-0.05 

50+0-02 
'J -0'-'- 0.02 

46.92+0.07 
65.5/58 
0.23 



2.033 
22.97:! 



3.333 
0.20"* 



+0.06 
0.06 
+0.11 



22.83_Q gg 

23.55«:- 

-,+0.16 
-0.12 
3+0.02 
-0.02 
3+0.12 
-0.08 
3 + 1.0 

46.93+6.08 



3.523 
0.10^ 
0.49^ 
2.3"* 



70.2/65 
0.32 

1 05+0.05 

22.8lto" 
23.08«-ig 

QQ fl'T"|~0.28 

20 

o /1O+0.20 
•^■■*°-0.26 
-,+0.01 
-0.01 
-+0.02 
-0.02 
3+0.06 
-0.06 
j+0.05 
-0.05 

46.91+0.08 



0.123 
0.15:! 
0.23:! 
0.48"* 



67.4/63 
0.33 



3+0.15 
-0.15 
„+0.12 
" -0.13 
00 on + 0.18 

ZO.OZ_g jg 

o trQ+0.25 
■J-'J='_0.17 
07+0-03 
'J-"' -0.06 
5+0.05 
-0.05 
1 4+0-6 
^■^-0.6 

47.10+0.16 



2.033 
22.76^ 



0.58"^ 



1.93 



56.1/65 
0.78 



2.00"^ 



3+0.10 
-'-0.10 

22.761°;}? 

23.50lo}| 

23.531°-}? 

3.4210;}? 

0.08l°;°^3 
so+0-05 

-+0.04 
-"-0.04 
1+0.03 
*-0.03 

47.08+0.09 



0.363 
0.54"" 



52.8/63 
0.82 



+0.07 
-0.05 
3+0.06 



1.91 



22 79^ 
23.421°;}? 

O cr, + 0.12 

■'•°^-0.25 
a+0.04 
-0.07 
3+0.06 
-0.07 
D+2.2 
'-1.0 

46.85+0.08 



0.283 
0.60j 
2.i 



89.6/60 
0.01 



+0.05 
0.05 
+0.09 
0.09 
+0.20 
0.22 



1.92 
22.79 
23.42 
< 23.5 

3 46+0-14 
■^•^0-0.13 

25+0-°2 
'J-^'^-0.02 
-+0.02 
-0.02 
a+0.05 
-0.05 
n =0+0.04 

46.87+0.07 



0.353 
0.48"* 



83.8/58 
0.15 



+0.09 
-0.07 
a+0.08 



22.78_Q Qg 
23.68l°;22° 
3 64+0-22 
3+0.03 
-0.06 
,+0.15 
-0.10 
-5+2.0 

46.95+0.10 



0.133 
0.62j 
2.6^ 



43.8/65 
0.98 



1+0.07 
-0.07 
a+0.06 
-0.08 
a+0.20 
-0.32 
3+0.20 
-0.36 
,+0.25 
-0.30 
1+0.04 
-0.03 
-,+0.08 
-0.08 
39+0-08 
"•■-'='-0.08 

0-581o;q7 
46.98+0.12 
41.3/63 
0.98 



1.913 
22.783 
23.48j 
23.40j 
3.52j 
O.li:* 
0.26j 



^ MODELl is a po-wcr-law with Galactic, neutral, and a one component outflow absorption (XSPEC model 
WABS*zwabs*mtable{windfull.tab}*PO'w); MODEL2 is a power-law with Galactic, neutral, and a two component outflow absorp- 
tion (xsPEC model wabs*zwabs*mtable{windvslow.tab}*mtable{windvfast.tab}*pow). 

'° The spectra fitted are the added (with HEASOFT ftools routine ADDSPEc) spectra of the FI chips (XISO, XIS2 and XIS3). For DBS 2 the 
XIS2 CCD was not operational, and only the XISO and XIS3 spectra were added for this observation. The BI spectra are taken with the 
XISl chip. 

^ Logarithm of the intrinsic (unabsorbcd) X-ray luminosity in the 2-10 keV rest-frame. These luminosities have been corrected for absorption 



(both Galactic and intrinsic). 

is referred to as slow and the one covering larger ener- 
gies (> 10 keV rest-frame) is referred to as fast. Typical 
models used to fit t he fast component contain a num 



ber of gaussian lines (phartas et al.]|2002 2007a 2009a 



Saez, Chartas fc Brandt ' 2U09p . We created tlie table 
model WINDVSLOW.TAB to ht the slow component and 
the table model WINDVFAST.tab to fit the fast compo- 
nent. Since table models windvslow.tab and windv- 
fast.tab are used to fit a relative small portion of the 
spectra we fix the velocity profile parameter to p = 1.0. 
Each component of the outflow has a velocity gradient 
At; with a minimum (wmin) and maximum velocity (t^max) 
that satisfy Wmax = ^^min + Aw. By fitting these table 
models to the spectra we constrain the minimum and 
maximum projected velocities of each component. Table 
model windvslow.tab is generated with Av having val- 
ues in the range 0.02-0. 42c (step intervals of 0.04c) and 
the parameters (wmin), log Nn and log ^ cover the same 
range and values as those used for producing table model 
windfull.tab. Table model windvfast.tab is generated 
with Aw having values in the range 0.04-0. 80c (step inter- 
vals of 0.04c for Aw < 0.24c and step intervals of 0.08c 
for Aw > 0.24c) and the parameters Wmax, log iVn and 
log ^ cover the same range and values as those used for 
producing table model windfull.tab. 



3.3.2. Results from fits to X-ray spectra of 
APM 08279+5255 

We fit the X-ray spectra of APM 08279-1-5255 using 
two different models. These two models assume a source 
spectrum consisting of a power-law attenuated by Galac- 
tic absorption (xspec model WABS). We a ssumed a 
Galac tic column density of 4.1 x 10^°cm~^ (Kalberla 



et al.||2005[ ). The first model (MODELl) contains an 
mtrmsic neutral absorber (xspec model ZWABS) to de- 
scribe the absorption in the low energy range of 1—4 keV 
rest-frame and an outflowing ionized absorber (xspec 
table model windfull.tab) to account for the X-ray 
BALs. In XSPEC notation MODELl is written as: 
wabs*zwabs*mtable{windfull.tab}*pow. The sec- 
ond model (M0DEL2) contains an intrinsic neutral ab- 
sorber (xspec model zwABS) to describe the absorption 
in the low energy range of 1—4 keV rest-frame and a 
two component outflowing ionized absorber (xSPEC ta- 
ble models windvslow.tab and windvfast.tab) to ac- 
count for the X-ray BALs. In xspec notation MODELl 
is written as: 

wabs*zwabs*mtable{ windvslow.tab}* 
mtable{windvfast.tab}*pow. 

We note that MODELl and M0DEL2 use a neutral 
absorber to describe the intrinsic attenuation in the low 
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TABLE 4 

Results from spectral fits of the Chandra and XMM-Newton observations of APM 08279+5255. 



Model" Parameter Values Epoch 1 Values Epoeh 2 Values Epoch 3 Values Epoch 4 Values Epoch 5 



MODELl. 



log {AfHi/cm-2) 
log {AfH2/cm-2) 
log $2 

p 

log L2-10 



, yt; + 0-04 
^•'^-0.03 
22 67+'^-°'' 
00 ^9+0.05 

ZO.OZ_g Qg 

-+0.08 
-0.12 
-+0.01 
-0.01 
,+0.03 
-0.04 
S+1.0 

46.83+0.06 
132.1/104 
0.03 



3.56: 
0.15^ 
0.43j 
3.8^ 



2 00+°"'^ 

r,g 79+0.04 
90 fi7+0.06 

ZO.D(_() gg 

-+0.09 
-0.10 
,+0.01 
-0.01 
1+0.07 
-0.06 
,+0.8 
-0.5 

47.01+0.07 
110.5/115 
0.60 



2.13^ 



3.35: 
0.13j 
0.64^ 
2.2"* 



fO.09 
-0.05 
99 yi +0.02 
zz. ( J-_n no 



-0.02 
-+0.08 
-0.05 
7+0.10 
-0.10 
a+0.02 
-0.02 
:i+0.03 
-0.03 
K+0-2 
-0.2 

47.15+0.08 
109.4/105 
0.36 



23.76: 
3.47j 
0.08j 
0.69j 
1.6^ 



9 OO+0.09 

Z-30_g Qg 

5+0.02 
'-0.02 
-+0.06 
-0.05 

r. + 0.05 

-0.05 
-+0.01 
-0.01 
,+0.03 
-0.02 
-+0.2 
-0.1 
47.48+0.12 

152.4/143 

0.28 



22. ( 
23.85;! 
3.36j 
0.05j 
0.73j 
1.5^ 



1 96+°°^ 
^•='"-0.03 
90 ri9 + 0.03 

ZO.UZ_g 03 

90 1 0+0. 09 
Z0.10_Q j^j^ 

o 90+0.12 
•^•^"-0.10 
,+0.03 
^-0.03 
,+0.02 



0.13"* 



31"* 
^•'-'^-0.02 
1 9+0.6 

46.91+0.04 
77.2/70 
0.26 



MODEL2., 



log (AfHi/cm-2) 
log (AfH2i/cm-2) 
log (A^H22/cm-2) 

log €2 
„(l) 



/C 

'max / 

(2) 



V 



log L2-10 



+0.02 
0.02 

+0.05 
0.05 
90 97+0.11 

ZCi.Z( _g^Qg 



1.77 
22.61 



a+0.12 
-0.12 
,+0.18 
-0.19 
-+0.01 
-0.01 
7+0.01 
-0.01 
,+0.02 
-0.02 
-,+0.01 
-0.01 

46.79+0.05 
116.8/102 
0.15 



23.08: 
3.32^ 
0.15^ 
0.17j 
0.32^ 
0.40"* 



9 04+0.05 

,+0.05 
^-0.05 

23.17l°;12 
23.47l«;l°3 

o 99+0.06 
'-'■^^-0.06 
13"^002 
'^■^•^-0.02 
1+0.04 
-0.04 



2.17 



+0.06 



22.63"* 



0.31"* 



< 0.12 
0.63"* 



,+0.04 
-0.03 

47.03+0.06 
107.5/113 
0.63 



22.67"* 



-0.07 
7+0.04 
-0.05 
23 60"'"°-'''' 

90 41+0.12 

q 07 + 0.08 

-0.08 

04+° °2 
*J-'J*-0.02 

0^40+0.05 



-0.05 
fO.lO 
-0.10 
,+0.03 
-0.03 

47.18+0.08 
106.5/103 
0.39 



0.43: 
0.69^ 



9 01 +0.08 
^•■-'J^-0.05 
,+0.03 
-0.03 

23.34lg-?^ 

,+o:o7 

-0.07 
90+0.04 
■J-^°-0.04 

06+° °^ 

'^•^^-0.05 
^+0.06 
-0.06 
,+0.02 
-0.02 

47.40+0.07 
152.5/141 
0.24 



22.63^ 



23.51^ 



0.24: 
0.73^ 



1 QQ+005 

j^-yy_o.o4 

1+0.04 
-0.04 

fJ(l-|-0.08 
ZO.ZZ_Q J^Q 

,+0.13 
-0.22 
q 9C-+0.15 
•-'•^''-0.11 

ii"^o oi 

^•^^-0.01 
,+0.02 
-0.02 
a+0.04 
-0.04 

'^•°"-0.03 

46.98+0.06 
66.2/68 
0.54 



23.01"* 



23.12"* 



0.31: 
0.58"* 



^ MODELl is a power-law with Galactic, neutral, and a one component outflow absorption (XSPEC model 
wabs*zwabs*mtable{windfull.tab}*pow); MODEL2 is a power-law with Galactic, neutral, and a two component outflow 
absorption (xsPEC model wabs*zwabs*mtable{windvslow.tab}*mtable{windvfast.tab}*pow). 



energy range of 1—4 keV rest-frame. In this work we also 
tried fits with an ionized absorber to describe the attenu- 
ation at 1—4 kcV; these fits did not result in a significant 
improvement over fits that used a neutral absorber. Ad- 



ditionally, as di s cussed in Saez, Chartas fc Brandt] ( 2009 ) ; 



Chartas et al.| ( |2009a[ ), the use of a neutral absorber 
to describe the absorption in the low energy range of 
1—4 keV rest-frame did not improve with the use of more 
complex absorbers even for spectral fits performed to 
deep X-ray observations of APM 08279+5255. We also 
note that given the high redshift of APM 08279+5255, 
our fits to the spectra of this source cannot adequately 
constrain the low-energy intrinsic absorption at 1—4 keV 
rest-frame. 

For M0DEL2 the current data cannot adequately con- 
strain the ionization parameters of both the slow and fast 
components of the outflow. We therefore set the ioniza- 
tion parameters of the slow and fast component to be 
equal in MODEL2.2" 

The results of the fits of these models to the Suzaku 
spectra of APM 08279+5255 are presented in Table [3] 
and to the Chandra and XMM-Newton spectra in Ta- 
ble H We note that in MODELl of Tables [3] and H 
log A'^Hi and log iVH2 are the best fitted column densities 



20 Notice that : 



L Saez, Chartas & Brandt (20091 we tried a sim- 
ilar model to MOUHjUI (model Ab'l'AK4, labie 5) and we found 
that the ionization state of the slow component could not be distin- 
guished for the ionization state of the fast component in the Suzaku 
observations. We have a similar case when we try MODEL2 on our 
eight X-ray observations without setting the ionization parameter 
of the slow and fast component as equal. 



of the intrinsic neutral absorber and the outflowing ion- 
ized absorber, respectively. In M0DEL2 of these same 
tables log A'^hi, log A^H21i log A^H22 are the best fitted 
column densities of the intrinsic neutral absorber, the 
slow and fast outflowing ionized absorber, respectively. 
In M0DEL2 of these tables the velocity superscripts (1) 
and (2) represent the best-fitted velocity parameters of 
the slow and fast component, respectively. 

In general the models that fit the high-energy out- 
flowing absorber (see Tables [s] and |4j) tend to have ion- 
ization parameters in the range JT2 < log ^ < 3.7. 
With the exception of the Chandra observations, there 
is no improvement in the flts when using M0DEL2 over 
MODELl. In Table [5] we show the statistical improve- 
ments based on the i^-test of fits using M0DEL2 over 
those using MODELl. For the Chandra observations 
the improvements are > 99.5% when we use M0DEL2 
over MODELl. There are also marginal improvements 
using M0DEL2 over MODELl for the Suzaku obser- 
vations especially in OBSl (see Table [s]). In general 
we will use M0DEL2 to estimate the outfiow param- 
eters because this model provides better constrains of 
the maximum velocity of the outfiow. However, in the 
case of the XMM-Newton observations we use MODELl 
to estimate the outflow parameters because these obser- 
vations show the widest X-ray BALs and provide bet- 
ter constraints of the velocity proflle parameter p. We 
mainly used MODEL 1 to fit the XMM-Newton spec- 
tra of APM 08279+5255 because only one broad ab- 
sorption trough is clearly observed at high energies. 
When we used Model 2 to fit the XMM-Newton spectra 
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TABLE 5 

Estimates of the improvement of fits to the spectra of 
APM 08279+5255 using MODEL2 over MODELl. 



OBSID Instrument F-statistic" null probability significance 



OBSl 


FI 


2.59 


4.8 


X 


10" 


2 


95.2% 


OBSl 


BI 


3.05 


1.5 


X 


10- 


1 


85.5% 


OBS2 


FI 


5.32 


2.4 


X 


10" 


1 


76.0% 


OBS2 


BI 


7.15 


1.4 


X 


10- 


1 


85.6% 


OBS3 


FI 


1.31 


2.8 


X 


10" 


1 


72.3% 


OBS3 


BI 


1.91 


1.6 


X 


10" 


1 


84.3% 


Epoch 1 


ACIS S3 


6.68 


1.8 


X 


10- 


3 


99.8% 


Epoch 2 


EPIC pn 


1.56 


2.2 


X 


10" 


1 


78.4 % 


Epoch 3 


EPIC pn 


1.40 


2.5 


X 


10- 


1 


74.9 % 


Epoch 4 


EPIC pn 














Epoch 5 


ACIS S3 


5.66 


5.32 


X 


10- 


-3 


99.5% 



The value on the left of the slash is the _F-statistic and is given by 



F = 



The value on the right of the slash represents 



the probability of exceeding the F-statistic based on the _F-test. 



of APM 08279+5255 the two outflow components were 
found to overlap in velocity space. This overlap occurs 



„(i) 



are less or similar to 



,(2) 



as pre- 



when values of Umax 
sented in Table IH 

Based on the spectral fits to the Suzaku observations 
(Table [3|, we find possible variability (> 2a) in the min- 
imum velocity of the outfiow between 0BS2 and 0BS3 

(wmin in MODELl and in MODEL 2) for both the 
BI and FI spectra. The evidence of this variability in the 
Suzaku spectra of APM 08279+5255 w as already ana- 



lyzed in Saez, Chartas & Brandt (2009| and in this work 
we conclude that its significance is at the 99.9% and 98% 
significance levels in the FI and BI spectra, respectively. 
This variability is also confirmed to exist between the col- 
umn density and velocity of the slow component when we 
compare epochs 0BS2 and 0BS3. 

We confirm the dramatic change in the maximum ve- 
locity of the fast com ponent of the Cha ndra observa- 
tions as shown in iChartas et al. (2009a). We find a 

(2) 

> 3(7 change of Wmax between Epoch 1 and Epoch 5 us- 
ing either MODELl or M0DEL2 (see Table E). We also 
confirm this change in velocity when we plot the con- 
fidence contours of the column density (log A'h22) ver- 
sus the maximum velocity of the fast outflow component 

f 2) I I 

) using M0DEL2. In Figure m we present confi- 
dence contours of log A'^h22 versus v\aix for Epoch 1 and 
Epoch 5. The significance of the change of log A^H22 ver- 
sus Vmax bctwecn Epochs 1 and 5 is > 99.9% with a null 
probability of < 10^'*. 

Using either MODELl or M0DEL2 we find a possible 
change in the maximum velocity of the outflow between 
the XMM-Newton observations of Epoch 2 and Epoch 
4, however, the significance of this change is only at the 
1— cr level (see Table |4]). We do find however significant 
changes in the minimum velocity and total column den- 
sity of the outfiow (log iVH2; Table |4]) between Epochs 2 
and 4. In Table [3] we show that this variability is signifi- 
cant at the > 2-cr level for fits using MODELl. In Fig- 
ure [7] we plot confidence contours of log A'h2 versus v„ 



(left panel) and log A'^h2 versus 



(right panel) for 



each XMM-Newton observation for fits using MODELl. 
We find that the significance of the variability between 



23.2 



E 
u 

\ 


23.0 


CM 
CM 
I 


22.8 


z 






22.6 






_o 


22.4 




Epoch 1 
Epochs 



0.4 0.5 0.6 0.7 
v(2) [cl 

max L J 



Fig. 6. — Contour plots of the column density (log (A''h22 /cm )) 
(2) 

versus maximum velocity (s^niax) of the fast component of the 
ionized outflow absorber. The contour plots have been calcu- 
lated for the Chandra observations Epoch 1 (continuous line) and 
Epoch 5 (dashed line) at the 68% and 99% level of significance 
using MODEL2 of Table [I] 



E 23. 



23.6 




en 
O 



23.4 



Epoch2 
■ Epoch3 
Epoch4 




0.05 0.10 0.15 0.5 0.6 0.7 0. 



Vmin [c] 



ax [c] 



Fig. 7. — Confidence contours of the total column density 
(log {Nii2/cni~^)) versus minimum velocity (i)min; left panel) and 
maximum velocity (fmax; right panel) of the high-energy absorber. 
The contour plots have been calculated for the XMM-Newton ob- 
servations Epoch 2 (solid line), Epoch 3 (dotted line) and Epoch 4 
(dashed line) at the 68% and 95% level of significance using 
MODELl of Table H 

Epoch 2 and Epoch 4 of log A^H2 versus Wmin is signifi- 
cant at the > 99.9% level and of log A'h2 versus Wmax is 
significant at the > 99% level. We also find variability in 
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vations of APM 8279+5255 (Fig ure 



1.8 2.0 2.2 2.4 

photon index (F) 



Fig. 8. — Maximum velocity (umax; upper panel) and to- 
tal column density of the outflow (log (A^jj/cm"^); lower panel) 
versus photon index (F) for our eight deep observations, timax, 
log (Afjj/cm"^) and F were derived from fits to the spectra of APM 
08279+5255 with a model that included an absorbed power-law 
and two component outflow (MODEL2 of Tables [s] and [ijl for the 
Suzaku and Chandra observations and a single component outflow 
for the XMM-Newton observations (MODELl of Table H. Errors 
shown are at the 68% confidence level. In each panel tne Suzaku 
data points are displayed as a single point that represents the 
weighted (by total counts) mean of the three Suzaku observations. 
Data shown with circles, a star, and squares are obtained from 
Chandra, Suzaku and XMM-Newton observations respectively. 

the minimum velocity of the outflow when we compare 
Epoch 2 and Epoch 3. This can be seen in Table g] and 
the left panel of Figure [Tj From the confidence contours 
(Figure m left panel) , we estimate the significance of the 
variability on log -/Vh2 versus Wmin between Epoch 2 and 
Epoch 3 to be > 99%. It is worth mentioning that the 
fits using M0DEL2 show an increase in the maximum 

velocity ~ 2— cr of the slow outflow component (wmlx) 
between Epoch 3 and Epoch 4. This change is consistent 
with the change found in the first component energy of 
the two gaussian absorption line model of Chartas et 
aLl (l2009a). Finally we note that the spectral fits to 



tEe xMM-Newton observations show values of p that are 
greater than one. Although we find indications of vari- 
ability of p, these changes are not significant given the 
poor constraints of this parameter. 

3.3.3. Mass outflow rate and kinetic energy of outflow of 
APM 08279+5255 

We use results from fits of our quasar outflow model to 
the spectra of APM 08279+5255 to constrain the mass 
outflow rate and kinetic energy injec ted by this out- 
flow into the s urro unding medium . In |Saez, Chartas fc 
Brandt] p009| and [Chartas et a l.| (|2009a| we also con- 
strained these properties of the outflow, however, these 
estimates were based on less realistic models (see [2.2]. 

A trend between Wmax and F based on the X-ray obser- 



upper panel) was 

recently found by [Chartas et al.| ( |20lI5a ). Such a trend 
is consistent with models that predict that quasar winds 
are driven by radiation pressure. Our fits with a more 
realistic quasar outflow model confirm the trend between 
'^max and F. We also find a possible trend between the 
total column density of the outflowing ionized absorber 
(log Nn) and F. As the lower panel of Figure |8] shows, 
we find that log Nh increases with F. The increase of 
the maximum velocity with F could be related to the 
fact that softer spectra have a s tronger radiative effec t 
on the wind (see 54] and §3.3 of Chartas et al. 2009a). 
On the other hana7a possible increase in colunin den- 
sity with F could be indicating that APM 08279+5255 
becomes a more effective radiative driving source as the 
spectrum of the ionizing source becomes softer. Addi- 
tional X-ray observations of APM 08279+5255 will show 
if these trends are significant. 

In order to calculate the mass outfiow rate (M) and 
the efficiency of the wind (ck) we have m odified the 
formulas used in Saez, Chartas fc Brandt ( [2009) and 



Chartas et al. 
ity gradient ot t 



2U09ap to include the modeled veloc 



corrections. 



le wind and include special relativistic 
The mass outfiow rate is given by (M = 

TVh 



4^/cmp|j^/o Hv(iVH)d7VH) 



M 



■' AR ^ 



+ 



1 



(4) 



and the wind efficiency, defined as the ratio of the rate 
of kinetic energy injected into the interstellar medium 
and IGM by the outflow to the quasars bolometric lumi- 
nosity, is (e/f = /(7 - l)rfM) 



In equations Q and ([5| fc is the global covering fac- 
tor. Nil is the column density, R is the radius, and AR 
is the thickness of the absorber. We note that for an 
absorber with constant non-relativistic speed {v/c ^ 1) 
equati ons ^ and ([s]) have the same form as the ones 



used in Saez, Chartas fc Brandt 



pogl) 

Saez, C 



and Chartas et al 



2009aD (see e.g., equation 4 of Saez, Chartas & Brandt 
^ The velocity profile assumed in our models is 
therefore equations (4| and ([5| 
fiov 



given by equation ([3|, therefore equations (4|) and 
contain the dependence of the estimated outflow proper- 
ties on the model parameters of the velocity profile. To 
obtain error bars for ex and M , we performed a Monte 
Carlo simulation, assuming a uniform distribution of the 
parameters fc, R, and R/ AR around the expected val- 
ues of these parameters, and a normal distribution for 



log A^H, 



and 



(described by the parameters in 



Tables |3] and [I]) . Specifically, we assume a global cover 
ing factor lying in the range fc = 0.1 — 0.3, based on the 



observed fraction of BAL quasars (e.g., Hewett & Foltz 
2003| [Gibson et ar] | 2009 [) and a fraction R/AR rang- 



ing from 1 to 10 based on current theoretical models of 

For a constant velocity, i.e. a p-type profile (equation 
[3| with p = and fmin = ''i then equation reduces to 
M = 4TTfc(R'^/AR)NnmpV. Also, since 7 1+v^/{2c^) (v/c < 1) 
equation K>\ becomes ex = Mv'^ /{2L]^al)- 



13 



TABLE 6 

Mass-outflow rates and efficiencies of outflows in APM 08279+5255 
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"•'J^-0.02 


1 2+1-1 


2+°'2 


r, 1+1.9 /r, r, + 2.1 

^■1-1.3/^-^-1.4 


0.2«-.?/0.3lg:^ 


Epoch2 
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3.9l^:|/4.3tl.? 
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^ Estimated values of the outflow properties of APM 08279+5255 based on a model that includes an absorbed power-law and a two 
component outflow (MODELS of Tables [s] and ^ for the Suzaku and Chandra observations and a single component outflow for the 



XMM-Newton observatio ns (M ODEL 1 of Table 



Mbh 

b 



10^ 



|4|. The values of M and e are obtained from equations l|4| and jsj 



assuming 



UJ 2 



(sec ^2.l\ and Lboi ^ 7 X lO"" /a^^Lq ( Irwin et al.|l998 Riechers et al.|2009 i. 

We provide two estimates of Af(tot) and eK(tot). The estimate on the left corresponds to the case where the angle between our 
line of sight and the outflow direction is zero while the one on the right corresponds to an angle of 10°. 

of MODELl for the XMM-Newton observations (see Ta- 
bles [s] and |4| . For the Suzaku and Chandra observations 
"we estimate M and for the slo-w and fast component 
of the outflow. The total mass outflow rate and efficiency 
is obtained by summing the contributi ons of t he slow and 
fast outflow component (M0DEL2, § 3.3.21 when mul- 
tiple components are required to fit the data. In the 
case of the XMM-Newton observations no summing is 
re quired since a one component wind model (MODELl, 
§ |3.3.2p provides better fits (in a statistical sense) to 
the spectra of APM 08279-^5255 (see § [3X2] ). In Ta- 
ble [6j assuming that the outflow is viewed along our 
line of sight, we show that during the observations of 
APM 08279-1-5255 the total mass outflow rate varied be- 
tween 2lt}^ - &btfi Mgyr-i (^l = 100) and the total 
efficiency varied between 0.2+"-^ — 2.3+^-g. However, the 

total values of M and are ~ 10% and ^ 50% higher 
respectively if we assume an outflow forming an angle 
of 10° with our line of sight (see Table |6|. In Table |6] 
we also show that an important fraction (> 10%) of the 
bolometric energy of APM 08279+5255 is injected into 
the surrounding galaxy through quasar winds. 

As Figure [8] shows the softer spectra probably result 
in a faster and more massive (higher total column den- 
sity) wind. This tendency can be seen more clearly in 
Figure [9] where we have plotted the total efficiency (^k) 
as a function of the photon index (F). We note that a 
small contribution to the error bars (< 20%) of tK arises 
from errors in the fitted parameters (i.e., Wmin, '^max and 
log A^h) while most of the error in ck is due to the un- 
certainty of _R, /c, R/ AR. 

The 0.5-8 keV flux of APM 08279+5255 varies by up 
to 50% in the observations of our study (see Table fll. 
These flux changes appear to be correlated to changes of 
the maximum outflow velocity (see Figure [lOl left panel). 
The observed changes in flux may be caused by changes 
in the luminosity of the ionizing source. By fitting our 
outflow model to the X-ray spectra of APM 08279+5255 



1.8 2.0 2.2 2.4 
photon index (F) 



Fig. 9. — Outflow eflicicncy (ex) of the outflowr versus photon 
index (F) for our eight observations of APM 08279+5255. tK and 
F were derived from flts to the spectra of APM 08279+5255 with 
a model that included an absorbedpower-law and two component 
outflow (MODEL2 of Tables [s] and |4| for the Suzaku and Chandra 
observations and a single component outflow for the XMM-Newton 
observations (MODELl of Table[4j. Errors shown are at the 68% 
confidence level. The Suzaku data points are displayed as a single 
point that represents the weighted (by total counts) mean of the 
three Suzaku observations. Data shown with circles, a star, and 
squares are obtained from Chandra, Suzaku and XMM-Newton 
observations respectively. 



quasar outflows ( Proga et al.||2000 2004). Based on our 
estimated maximum velocities (fmax ~ 0.7c) (see e.g., 
Chartas et al.||200"9a ) and the fast variability of the out- 
flow we expect that R will be close to the Schwarzschild 
radius (Rs)- Therefore in the Monte Carlo simulation 
we allow R to vary between 3i?s and 15i?s. 

In Table [e] we show the mass outflow rate (M) and the 
efficiency (ex) of the wind in each observation. The pa- 
rameters used in equations ^ and ([H) to derive M and 
eK were obtained from spectral fits of M0DEL2 to the 
Suzaku and Chandra observations, and through the use 
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'0.5-8 



2-10 



Fig. 10. — Maximum velocity (tJmax; upper panel) versus ab- 
sorbed flux in the 0.5—8 keV band (/o.5— left panel) and in- 
trinsic luminosity in the 2—10 keV rest-frame band (L2— lo; right 
panel) for our eight deep observations, timax and L2-10 were de- 
rived from fits to the spectra of APM 08279-1-5255 with a model 
that included an absorbed power-law and two component outflow 
(MODEL2 of Tables Js] and ^ for the Suzaku and Chandra ob- 
servations and a single component outflow for the XMM-Newton 
observations (MODELl of Table H. The absorbed fluxes in the 
0.5—8 keV band are from Table Errors shown are at the 68% 
confidence level. In each panel the Suzaku data points are dis- 
played as a single point that represents the weighted (by total 
counts) mean of the three Suzaku observations. Data shown with 
circles, a star, and squares are obtained from Chandra, Suzaku and 
XMM-Newton observations respectively. 



we obtain the intrinsic (unabsorbed) 2—10 keV rest- 
frame X-ray luminosity of the ionization source (^2-10) 
(see Tables|3]and|4|. We find that the maximum velocity 
of the outflow is possibly correlated to t he i ntrinsic lumi- 
nosity of the ionizing source (see Figure 10 right panel). 



From the radiative outflow wind equation we obtain that 
the terminal velocity of the wind is Voo oc ^^TfL, where 
L is the luminosity producing the radiative driving and 
Fy is the force multiplier (see §4 for an analysis describ- 
ing the dependence of F^ with the photon index). This 
equation predicts that an increase in the luminosity of 
the ionizing source should result in an increase of the 
terminal velocity of the outflow as observed. The Suzaku 
observations of APM 08279-1-5255 do not show any sig- 
nificant changes of F, Wmax, the total colum n de nsity, and 
(see Table [S]). In Figures [S) [9] and 



10 we there- 



L2-10 

fore chose to represent the Suzaku observations with one 
data point. The Suzaku data point in each figure is the 
weighted mean of the three Suzaku observations. 

As a final comment, the ionization parameter obtained 
from the spectral fits is based on the assumption that the 
incident spectrum is a pure power-law with a fixed slope 
(a = — 1). In reality the spectra are likely more complex 
and variable than assumed in our model. Our simple ap- 
proach provides reliable constrains on the column density 
and velocity of the outflow. However, since the ionization 
parameter depends on the SED from 1—1000 Ry (see Ap- 
pendix|A]) we do not expect that fits of our quasar outflow 
model to the spectra of APM 08279 -1-5255 to constrain 
the ionization parameter since our model assumes a fixed 
SED. What our model does constrain is the temperature 
of the absorber which is a reliable indicator of the degree 
of ionization of the outflowing gas. 

For example, for a fixed column density and velocity profile 
we have checked that approximately the same absorption features 
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Fig. 11. — Modi fied Mathcws-Ferland SEDs used to calculate 
force-multipliers in §4.1| The x-axis is the logarithm of frequency 
(in Ry units) and the y-axis is in units of the logarithm of the fiux 
times the frequency (arbitrary units). The left and right panels 
correspond to two modified versions of the Mathews-Ferland SED. 
In the left and right panels F = 2.1 and F = 1.7 with cvox = —1.52. 
The middle panel corresponds to the standard Mathews-Ferland 
SED, i.e. F = 1.7 with ctox = —1.41. In each panel the dotted line 
is the unabsorbed SED and the solid line is the shielded SED (with 
log C,h~2.9). 

4. INFLUENCE OF THE SED ON THE DYNAMICS OF THE 
OUTFLOW 

In section 3 we presented a quasar outflow model that 
was fit to X-ray observations of APM 08279-1-5255 in or- 
der to constrain the column densities of the absorbers, 
their ionization state and their outfiow velocities. The 
outflow models used to fit these X-ray data assumed 
that the incident SED is a power law with a spectral 
index of a = — 1. In this section we attempt to under- 
stand the influence of the incident SED on the dynam- 
ics of the outflow. We vary the SED to obtain insight 
into the dependence of the force multiplier on the inci- 
dent SED, however, we do not fit any data with variable 
SEDs. This section extends t he analysis presented in 
§3.3 from Chartas et al. (2009a I by providing an analysis 
of force multipliers as a function of the sp ectral changes 
of a Mathews-Ferland SED, while, in the [Chartas et ah] 
(j2009a) paper we performed the same analysis but using 
a pure power-law SED. Here we focus on the new results 
that are based on more realistic SEDs. 

4.1. The dependence of the force multiplier on the 
photon index and aox for a Mathews-Ferland SED. 

As suggested in Chartas et al. 2002, 2009a and Saez et 
al. 2009 it is plausible that the driving force on the high- 
energy absorbers is produced initially by X-rays. The 
short term variability time-scale of the X-ray BALs of 
APM 08279-^5255 suggests a launching radius of < 10i?s 
and recent studies of AGN employing the microlensing 
technique indicate that the X-ray emission region of the 
hot corona in AGNs is compact with a half-light radius 
of a few Rs and t heir UV regions are roughly a facto r 



of ten larger (e.g., Morgan et al. 2008 Chartas 



2009b) 



Therefore UV radiation is not expected to contribute ini- 
tially at small radii to driving the X-ray absorbing out- 
are generated if we use an incident power-law SED with a = —1.0 
and log ^ 3.4, an incident power-law with a = —0.6 and log ^ 
2.7, an incident power-law with a = —1.2 and log § ft: 3.8, or an 
incident Mathews-Ferland SED with log ^ 4.7. We note that 
the average temperatures of the gases producing these absorption 
profiles are approximately independent of the incident SEDs and 
lie within the range of 6.6 < log T < 6.8. 
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Fig. 12. — The bound-free, Mq, and bound-bound, A/l, com- 
ponents of the force multiplier are shown as a function of the ion- 
ization parameter. Force multipliers are calculated for Mathcws- 
Ferland SEDs with photon indices of F = 2.1 (solid lines) and F = 
1.7 (dotted lines). In the left panel we have assumed no absorbing 
shield, whereas, in the right panel the soft and hard SEDs have been 
attenuated by a warm-absorber shield with log (Afjj/cm"'^) = 23 



and log gj. 



2.9 on its illuminated face. The SEDs used are 



described by the middle and right panels of Figure |11| 

flow. However, as the outflowing absorber gets further 
away from the source the contribution of the UV pho- 
tons to the driving force wiU increase relative to that of 
the X-ray photons. To investigate the driving mecha- 
nism of the wind we calculated the force multipliers for 
SEDs that extend to radio wavelengths and estimated 
the effect of changes in F. In Chartas et al. (2009a) we 
described the dependence of the force multiplier on the 
photon index assuming the SED is a power-law extend- 
ing from the UV (or 1 Ryd ~13.6 eV) to hard X-rays 
(or lO"' Ryd ^100 keV). In those calculations we assume 
SEDs with two different values of the power-law pho- 
ton index. A "soft" (E = 2.1) and a "hard" (F = 1.7) 
SED. In this work we extend our simulations to include a 
standard and a slig htly modified version of the Mathews- 
Ferland SED (Math ews fc Ferland||1987l ). For each SED 
we calculated the continuum {Mq) and the line (Ml) 
components of the force multiplier. depends on an 
additional parameter, t, which is commonly referred 
to as the "effective electron optical depth" and encodes 
the dynamical information of the wind in the radiative 
acceleration calculation (see Appendix |e]) . For our cal- 
culations we have assumed log t — —7. 

The standard Mathews-Ferland SED is characterized 
by ttox = —1.41 and by a power-law spectral index 
of ax = -0.7 (F = 1.7) which extends from 27 Ry 
(-0.36 keV) to 7.4x10^ Ry (-100 keV). The modified 
Mathews-Ferland SED is similar to the standard one but 
with ax — —1-1 and aox = —1-52 (see Figure [ll] to see 
the different SEDs used). In the left panel of Figure 12 
we show that for low ionization parameters (log ^ 
the UV dominates the driving of the wind, therefore in 
this regime there is no major change in the force multipli- 
ers when we compare the soft and hard Mathews-Ferland 
SEDs. However at log £, }^2, the line force multipliers for 
the soft Mathews-Ferland SED case are larger than those 
for the hard SED case. This effect is produced because 
in the case of a Mathews-Ferland SED with a soft X-ray 
spectrum (F = 2.1) the outflowing gas does not become 
very ionized at high levels of the ionization parameter 

The dimensionless optical depth is t = n^aTV-n^/{dv /dr), 
where, ne is the electron number density, cx is the Thomson cross 
section and Vtjj is the thermal velocity of the gas. The line force 
multiplier increases with decreasing t. 



Fig. 13. — The bound-free, Mq, and bound-bound, A/l, com- 
ponents of the force multiplier are shown as a function of the ion- 
ization parameter. Force multipliers are calculated for a Mathews- 
Ferland SED with Qox = —1.52 (solid lines) and Oox = —1.41 
(dotted lines); in both SEDs F = 1.7. In the left panel we have 
assumed no absorbing shield, whereas, in the right panel the soft 
and hard SEDs have been attenuated by a warm-absorber shield 
with log {Nn/ cm~'^) = 23 and log ^^j^ = 2.9 on its illuminated 
face. The SEDs used are described by the middle and left panel of 
Figure 

compared to the case of a Mathews-Ferland SED with 
a hard X-ray spectrum (F = 1.7). We conclude that an 
ionizing source that is relatively soft in the X-ray band 
will result in a force multiplier that remains relatively 
large even at ionization levels of the outflowing absorber 
of log C ^ 3. This prediction is consistent with the ob- 
served correlation of maximum outflow velocity with the 
X-ray photon index. We modified the SEDs by including 
a warm-absorber shield with log (iVn/cm"^) = 23 and 
log ^gj^ = 3.1 (Figure 12 right panel; see also Figure [TT| 
to see the shielded SED). The motivation of the inclu^ 
sion of this absorbing shield is based on the fact that 
both observations and theoretical models of radiatively 



driven winds suggest it s existence (e.g., Murray et al. 



1995||Proga et al ||2004[ ). We note that the observed ab 
sorption at rest-frame energies < 2 keV appears to be 
produced by gas of low ionization parameter and could 
therefore be associated with this proposed shield. The 
main effect of including a shield is a larger increase of 
the line force multipliers for the soft SED case (F = 2.1) 
compared to the increase produced by the hard SED case 
(F — 1.7) for log ^ > 2. We also find that by including a 
shield we obtain significant line force multipliers at ion- 
ization levels that are larger than the ones possible with 
no shield present. This effect is especially important in 
the soft SED case. 

In the process of changing the X-ray photon index from 
the default value of F = 1.7 (hard-spectrum) to F = 2.1 
(soft spectrum) in the Mathews-Ferland SED, we also 
changed aox from —1.41 to —1.52. To test the sensi- 
tivity of the force multiplier to changes of the aox of the 
SED, we estimated the force multiplier for the case where 
aox changed while the X-ray spectral slope remained con- 
stant. In order to produce this effect we increased the 
blue bump of the Mathews-Ferland SED by changing the 
default value of aox from —1.41 to —1.52 (to see the SEDs 
used see Figure 11). As shown in the left panel of Fig- 
ure [13] the force multiplier is more sensitive to aox than 
F. This difference in sensitivity is associated with the 
fact that the UV covers the range where most of the 
absorption lines lie. We also modified the SEDs by in- 
cluding an absorbing shield with log (iVn/cm^^) = 23 
and log ^gjj = 3.1 (see Figure 11). In the right panel of 
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Figure 13 we show that the effect of adding this warm- 
absorber shield is an increase of the force multipher for 
soft SEDs even at high ionization levels. 

We note that we used a moderate amount of shielding 
for the estimates shown in the left panels of Figures [12] 
and 1131 The combination of moderate attenuation and 
high force multipliers in this type of shield is expected 
to pr oduce high velocity outflows (Chelouche & Netzer 
2003 1 . If we use an absorbing shield with a relatively low 



ionization level we find that the large absorption from the 
shield will attenuate a large fraction of the UV and X-ray 
photons, and therefore prevent the outflow from reaching 
high terminal velocities. 

In general, the presence of an absorbing shield will in- 
crease the differences between the line force multipliers 
for the soft and hard SED cases. These differences will 
be important for shields with ionization parameters in 
the range of 1.0 < log ^gj, < 1.8 for the power-law SED 
case and will be important for shields with ionization 
parameters in the range of 1.2 < log < 3.1 for the 
Mathews- Ferland SED case (see Appendix IfI). For val- 
ues of the ionization parameter of the shiela over these 
boundaries, the shield will be transparent to incident ra- 
diation and the force multipliers will be si milar to one s 
with no shield present (left panels of Figures 12 and 13). 
For values of log < 1.0 the large absorptio n produced 
by the shield wi ll result in a weak outflow (see, Chelouche 
fc Netzer|[2003| . 

If the F— Umax trend is confirmed with additional obser- 
vations we plan to produce a more sophisticated model 
that will include more realistic kinematic, ionization and 
absorption properties of the outflow. If the F — u,nax re- 
lation is real we would also like to know if the changes 
in F are accompanied by changes of aox- Studies of 
correlations between outflow properties and the SED of 
the ionizing source will provide crucial information re- 
garding the driving mechanism of the outflow. Simul- 
taneous multiwavelength observations of sources such as 
APM 08279-^5255 that show fast outflows will therefore 
be essential for the study of quasar outflows. We note 
that for the SEDs analyzed in this section the contri- 
bution to the force multiplier from highly ionized iron 
lines such as Fe xxv Ka (Is^ — ls2p; 6.70 keV) and/or 
Fe XXVI (Is — 2p; 6.97 keV) is relatively smaU. How- 
ever, we observe significant blueshifts of the X-ray BALs 
of APM 08279-1-5255 indicating highly ionized material 
moving at ncar-relativistic velocities. One possible inter- 
pretation is that we are observing absorption from out- 
flowing highly ionized material that was accelerated at a 
previous time when it was at a lower ionization level. We 
note that a better understanding of quasar outflows will 
have to await for the development of more sophisticated 
numerical codes that combine MHD simulations of the 
outflowing material and the photoionization properties 
of the wind. Such an analysis is beyond the scope of this 
paper. 

5. CONCLUSIONS 

We have re-analyzed 8 long exposure X-ray observa- 
tions of APM 08279-^5255 (2 Chandra, 3 XMM-Newton, 
and 3 Suzaku) using a new quasar-outflow model. This 
model is based on cloudy simulations of a near- 
relativistic quasar outflow that were used to generate 
XSPEC table models. 



The main conclusions from this re-analysis are: 1) In 
each observation we have found X-ray BALs that have 
been detected at a high level of significance. Specifi- 
cally, we find strong and broad absorption at rest-frame 
energies of 1 — 4 keV (low-energy) and 7—18 keV (high- 
energy). 2) We confirm that the medium producing the 
low-energy absorption is a nearly neutral absorber with 
a column density of log (iVH/cm~^) ~ 23. The medium 
producing the high-energy absorption appears to be out- 
flowing from the central source at relativistic velocities 
(between — 0.7c) and with a range of ionization param- 
eters (3 < log ^ < 4). The maximum detected projected 
outflow velocity of ~ 0.7c constrains the angle between 
our line of sight and the wind direction to be < 22° . The 
short time-scale variability (< week in the rest-frame) of 
the high-energy absorption implies that X-ray absorbers 
are launched from distances of ^ 10i?s from the cen- 
tral source (where Rs is the Schwarzschild radius). 3) 
We also confirm trends between the maximum projected 
outflow velocity (umax) with the photon index (F) and 
the intrinsic (unabsorbed) X-ray luminosity at the 2—10 
rest-frame band; i.e. softer and more luminous X-ray 
spectra generate larger outflow terminal velocities. We 
also find that the total column density (-/Vh) of the out- 
flow increases with F. The trends found suggest that the 
wind becomes more powerful as the incident spectrum 
becomes softer. This tendency is confirmed through es- 
timates of the efficiency of the wind. Additionally, we 
estimate that a significant fraction (>10%) of the total 
bolometric energy over the quasar's lifetime is injected 
into the intergalactic medium of APM 082794-5255 in 
the form of kinetic energy. 

In this work we modeled the spectrum of the cen- 
tral ionizing source with a power-law SED and with 
Mathews-Ferland SEDs and found that variations of the 
X-ray and UV parts of the SEDs will produce impor- 
tant changes in the strength of the radiative driving 
force. These results support the observed trend found 
between the outflow velocity and X-ray photon index 
in APM 08279-1-5255. In general we find as expected 
that the presence of a mode rate ab sorbing s hield result s 
in more powerful outflows (Chelouche & Netzer 2003). 
Specifically, we find that shields with column densities 
of log (iVn/cm"^) ^ 23, covering an optically thin out- 
flow, provide a significant increase in the driving force 
when their ionization parameters are in the range of 
1.0 < log fsij < 1.8 for the case of a power-law SED 
and in the range of 1.2 < log ^ 3.1 for the case of 
a Mathews-Ferland SED. Therefore the strength of the 
radiative driving force depends critically on both the col- 
umn density of the shield and its ionization level. A con- 
firmation of the results found in our simulations of quasar 
outflow will require new deeper X-ray multi-wavelegth 
observations of quasars that contain clear signs of fast 
outflow. Such observations will allow us to correlate the 
properties of the outflow with properties of the SED and 
thus test our predictions. 
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Fig. A.l. — Limiting column density at which a layer will be considered thin. In this case we arbitrarily define a layer as thin when the 
ionization parameter decreases after crossing the layer less than 0.05 dex due to absorption. The solid and dotted lines represent a soft 
spectrum power-law (F = 2.2) and hard power-law spectrum (F = 1.5), respectively. Both spectra are defined between 1 Ry and lO** Ry 



APPENDIX 
IONIZATION PARAMETER. 

There are two definitions of the ionization parameter that are broadly used in the scientific literature. The first 



definition, which will be mostly adopted in this work, is the ionization parameter of Tarter, Tucker & Salpeter ( 1969 1 
given by 

e=^ = — / F,d,^= (Al) 

where Lj is the ionizing luminosity, Fj is the ionizing flux, Fi, is th e incident flux, u h is the hydrogen density, and 
r is the source-cloud separation. A second definition of the ionization ( Davidson|1977 ) parameter is given by the ratio 
of photons that can ionize hydrogen to the number of hydrogen atoms in a spherical layer at a distance r from an 
illuminating point source, i.e., 

U= ^ , where Q= -^dv. (A2) 

Airr-^cnH Jmy 

We note that for a pure power law (i.e., F^ cx v°'), the analytic expressions to convert between ^ and U (valid for 
a < 0) are, 

f log f/ 4- 1.754 a = -1 

log ^ = j log U + 0.914 log (-(1000"+i - 1)^) ^^^^ 



WHEN CAN AN ABSORPTION SLAB BE CONSIDERED THIN? 

The fraction of flux absorbed in the rest-frame of a layer of gas depends on the ionization state of the absorber, 
its column density, and the SED of the illuminating source. Our quasar outflow code calculates the absorption 
profiles resulting from absorption of the central source spectrum by ionized gas in an accelerated outflow. One of 
the assumptions of our quasar outflow model is that the absorption layer producing the observed X-ray BALs has a 
plane-parallel geometry with a thickness that is much smaller than the distance from the central source^''. We therefore 
expect any decrease in the ionization parameter across an absorbing layer to be predominantly due to the absorption 
of the flux in the layer rather than to the 1/R^ decrease of the ionization parameter. 

We will refer to a layer as optically thin if the change between the ionization parameter on the illuminated side and 
the dark si de of the layer is less than 0.05 dex.^^ 

In Figure we show the column density of a thin layer that will produce a decrease of 0.05 dex of the ionization 
parameter across this layer as a function of the ionization parameter at the front surface of the layer. We plot this 
relation for power-law SEDs with a photon index of F = 1.5 (dotted line) and with F = 2.2 (solid line). The results 

•^^ Let's assume a relatively faint AGN with Lj ^ lO^^erg s^-'^ with an absorbing layer of log (A'jj/cm^^) = 22, njj = lO^cm"'^, 
and an ionization parameter of log § = 3 on the illuminated side of the layer. Since ng ^^'^ cm; therefore 

AR/R ~ 0.3 rig ^^"^ . Since this is an extreme case (faint AGN and highly ionized layer), in general we expect AR/R <C 0.1. 

If we define and if as the ionization parameters at the illuminated front and back sides of the layer, respectively, then our definition 
of a thin layer can be written as log (C//Ci) < —0.05. In addition if we assume that the incident flux on layer is F^, we have that 
log = log(/^p Fve~'^" dv/ f^^ F^du) ~ log(l — (r)) ~ — (r) /InlO, where = I Ry, vi = 1000 Ry and (r) is the average opacity 

through the layer between uq and i/i. Therefore, our definition of a thin layer is equivalent to a layer with an average opacity of (r) < 0.1. 
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Fig. a. 2. — 0.3-13 keV spectra of the incident (full-line), transmitted (dashed-linc) and emitted (dotted-line) continua for different 
ionization levels of the outflowing material simulated with our quasar outflow code. In Panels 1-3 (Top, left to right) the ionization levels of 
the outflowing material are log 2.75, 3.0 & 3.25, respectively. In Panels 4-6 (Bottom, left to right) the ionization levels of the outflowing 
material are log 2.75, 3.0 & 3.25, respectively. For all these simulations the total number of layers is Ni = 100, the total column density 
of the outflowing material is A'^h = 3 x 10^^ cm"'^, the initial gas density is lO^'^ c m~^ and the velocity of the outflow is assumed to vary 
linearly between and 0.3 c. These results are in good agreement with those of |Schurch fc Done| | [2007| l as evident from a comparison 
between this flgure and their Figure 7. 

show that the upper Umit of the column density of an absorbing layer required for this layer to be considered thin, 
increases with ionization parameter. Soft spectra have lower thresholds than hard spectra, especially for low values of 
the ionization parameter (log ^ < 2). However, for highly ionized media the threshold of the soft and hard spectra are 
similar, in particular for (log ^ > 3) the limiting column density is log (A^n/cm^^) ^ 23. 

NEAR-RELATIVISTIC QUASAR OUTFLOW CODE: ASSUMING AN OUTFLOW WITH A CONSTANT AV BETWEEN 

LAYERS 

Since our library of cloudy simulations was generated for a specific turbulent velocity we adjusted the relative speed 
between layers i and « -I- 1 to be constant (i.e. lS.v[ = (I'i+i — Wi)/(1 — ViVi+i/c^) constant). Assuming Ni layers, we 



initialize the parameters defining our profile by: 1) Deriving the velocity profile, which satisfies vq 



Vn 



(Ati^ = Aw'). 2)Deriving iVn^ (i=0,..,A^;) through the use of A^nf+i — -^nf = ^Vi/k; this result comes from v{Nii) 
Wmin + kNu^ {k = (wmax — Vmin)/]^HT and Nht is the total column density of the outfiow). Using the last expression 

_i starting from A'^ho = 0- 3) Calculating the average velocity Vmi = {vivi+i) + v{ri))/2, 



A^Hi, the density rii = n{vmi), (i=0,..,A^i_i) and the radial steps Ar^ = ANn/m 



we calculate Nm from 
the column density AA^Hi = 
of each layer. 

We note that our assumptions are valid as long as Aw < Wtuib\/l2 (see 53.2.3). In our simulations we chose Avi ~ Wturb 
and Wturb = 1500 km s^^ ^ 0.005c . In the case of Vmin — 0.1c and Wmax = 0.7c we would have Ni ^ 100 layers. 



NEAR-RELATIVISTIC QUASAR OUTFLOW CODE: COMPARISON WITH SCHURCH & DONE (2007) 



In this appendix we compare our quasar outflow code to a similar code presented in Schurch & Done (2007). The 
goal is to estimate the level of agreement between two independent quasar outflow codes and at the same time evaluate 
the limitations of each approach. We first generated absorption profiles for a linearly accelerated outflow {p — 1) with 

10^^ cm-3, and values of log ^ = 2.75, 3, 3.25, 3.5, 3.75 and 4. In 



0, 



0.3c, A^H = 3 X 10^3 cm 



no 



order to make a fair comparison with the results of Schurch & Done ( 2007 ) we ran a large number of simulations using 
the photoionization code xstar and assumed a power-law model with a = —1.4. 

In each layer we subtract the absorption an add the emission. Since XSTAR does not include a "Compton scattered" 
component for the emission, we add this component in a similar way as the one described in Schurch & Done 2007. 
The Compton scattered emission component is obtained by calculating the Thomson sc attered flux of each lay er^^ 



Schurch fc Done| p007| ) the 
Dal covering factors equal to 



where the incident spectrum is the output from the previous layer. Since in the case of 
calculations were made assuming spherical symmetry of the multilayer absorber, with glo 
one, we use similar assumptions to obtain the emission of each layer. We calculated the emission from a spherical shell 
by adding the spectral contribution from each solid angle of the layer. The emission from each solid angle of the shell 
is corrected for Doppler beaming of material outflowing at a projected velocity with respect to a fixed line of sight. 

The spectra of the incident, absorbed and emitted continua for each of the initial values of the ionization parameters 
used in these calculations of absorption line profiles are presented in the six panels of Figure 
similarity of Figure |A.2| with Figure 7 of Schurch & Done 2007. The small differences 



A. 2 We find a remarkable 
"ound in our results when 



^® The Thomson scattered flux is the incident flux multiplied by 1 
Thomson cross section. 



AJVhctt where AAfjj is the column density of the layer or is the 
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Fig. a. 3. — Emitted (upper panel) and absorbed (lower panel) spectra in the enerffv range where we detect iron absorption (~7-9 keV) 
for three different values of the exponent of the p-type velocity profile, (see equationls] of p = 5.0, 1.0 and 0.2. For these comparison 



tests we set v„ 



0, Vn 



: 0.3c, iVn = 3 X 10^3 cm" 



no 



10 cm 3, and log ^ = 3. We find agreement between our code and that 



of Schurch & Done 2007 as shown in the comparison between the absorption and emission profiles presented in Figure [A. 3| and Figure 4 of 
[Schurch fc Done| | |2007l l . 



compared to those of Schurch & Done (2007) are mainly associated with differences in the details of the calculations. 
Specifically, Schurch & Done 2U07 pass the output flux from one layer to the other internally^^, while in our approach 
we perform separate photoionization runs from one layer to the next. Possibly an important factor that contributes 
to the discrepancies between the two codes is that in our calculations we approxim ate the incident flux on each layer 
as being the unabsorbed source SED. On the other hand, Schurch & Done (20071 use the absorbed SED to perform 
the photoionization calculations in each layer. We expect our approximation of using the unabsorbed source SED as 
incident on all layers to break down at low ionization levels where absorption th rough a layer may signifi cantly change 
the SED. This explains why our results begin to noticeably differ from those of [Schurch fc Done] (2007) at log C ^ 3. 

In Figure A. 3 we plot the absorbed and emitted spectra in the energy range where we detect iron absorption (^7- 
9 keV) for three different values of the exponent of the p-type velocity profile, (see equation Isl ^ of p = 5.0, 1.0 and 



0.2. For these comparison tests we set Vn 



0, wniax = 0.3c, iVn = 3 X 10^^ cm-^ no = lO^-^ cni^^ and log ^ = 3. We 



find agreement between our code and that of S church fc Done 2 007 as shown in the com parison between the absorption 
and emission profiles presented in Figure A. 3 and Figure 4 of Schurch fc Done (2007). Our emission spectra for the 
ca se of a velocity p r ofile w ith an exponent of p = 5 differs from those derived by Schurch fc Done 2007 (see Figure 4 
of Schurch fc Done ( 2007[) ) mainly because in their simulations the turbulent velocity is set to be variable in order to 
model fine structures in the spectra. In our simulations we fix the turbulent velocity in the outflow, and therefore we 
are not able to reproduce flne details (< 1500 km s^^) in the spectra. This is not a significant problem since we are 
mainly interested in simulating the broad absorption lines detected in the X-ray spectra of BAL quasars.^* 

We conclude that our quasar outflow code provides results that are very similar to those found in the independent 
code of Schurch fc Done 2007 as long as the SED does not change significantly across the multilayer. For our study 
we have mostly selected quasar outflows of highly ionized material (i.e., log ^ > 3.0 ) and therefore expect the SED 
of the central source not to be significantly attenuated across the outflowing multilayer. As a flnal remark, we note 
that in our quasar outflow code the photoionization calculations at each layer use a library of preexisting runs that 
can be performed either with cloudy xstar or any other photoionization code. This presents an advantage over the 
quasar outflow code XSCORT in the sense that we can compare photoionization calculations performed with different 
photoionization codes which might incorporate slightly different atomic physics. Additionally, the major advantage of 
using our outflow code over XSCORT, is a significant shortening in the computing time required for the simulation of 
an absorption profile. 



27 In | |Schurch fc Done|2007| l the authors use the same photon field XSTAR propagates internally, therefore their approach is consistent 
with a single, constant density, XSTAR run, in the absence of an external velocity field and constant gas density throughout the chain. 

2* A typical energy resolution of a X-ray telescope (at energies between 1—10 keV) is ~ 100 keV, therefore for a Fe xxv Ka (Is^ — ls2p; 
6.70 keV) line we will not be able to resolve features with Doppler broadenings < 4000 km s~^. 
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FORCE MULTIPLIER CALCULATIONS 

In this append ix we derive a formula for the force multipl ier and compare our expression to ones reported in earlier 
studies (see e.g., |Arav et al.||1994[ |Mihalas fc Mihalas||1984[ ). 

Force multipUer calculation 
For a medium that absorbs radiation, the force per unit volume on the medium is 

fv = \j aMv, (El) 

where a^{cm^^) is the absorption coefficient and Fi, is the incident flux on the medium. Assuming that the flux 
originates from a point source with luminosity at ?' = the acceleration is given by 



dv{r) 1 
dt 4:Trr^p{r)c 



a^L^e-^-'du. (E2) 



The factor e corresponds to the attenuation factor at radius r, [T^{r) — a^{r')dr'). Assuming that the central 
source has a mass of Mg the force equation can be re-written as: 

— ^ ^ , ''°' Mir) ^ with Mir)^^ , \ , , E3 

where Lboi is the bolometric luminosity, ut is the Thomson cross section and M(r) is the force multiplier. M{r) 
can be interpreted as the ratio of energy absorbed by the medium to the energy absorbed by the electrons. 

As a simple example we assume that the medium is a spherical shell with inner and outer radii of R and R + Ai?, 
respectively. For this case {M{r)) = J M{r)dr / AR; and therefore, the wind equation (E3) is approximated as 

dt 1.2TOp47rr2c ^ r2 ^ ^ 

In this last stage we have assumed a high ionization state for the gas at solar metallicities and t here fore pr/ne{r) ~ 
1.2mp, where m-p is the mass of the proton. Assuming that i'(-Rmin) = the solution of equation (E4) is 

where is the Eddington Luminosity (i.e. LEdd = i'^cGiimpM^ y 



To calculate the force multiplier of a line we start from equation ( |E3[ ). Let's assume that the absorption line has a 
constant opacity, frequency vi and a width Az^th = vivtii/c (where wth is the thermal speed of the absorbing atoms^^). 
The characteristic distances of radiative interactions are of the order of I ^ Wth/(Vi;). Hence for radially streaming 
radiation in an expanding medium the optical thickness is, 

n = aivth/{dv/dr) (E6) 

where a; is the absorption coefficient of the line. We define rji ~ ^"^^ and introduce an equivalent electron optical 
depth scale 

t = Ti/rji = neO-TVth/idv/dr). (E7) 
Therefore the mean force multiplier is (A/(r)) = J M{r)dr/l, i.e., 

(Mir)) . = (E8) 

Summing up over all absorption lines we obtain 

M{t) = - — > 7 w - — > L^^Avth,i mm 77,, - . (E9) 

-t^bol J t Lho\ ^ \ t J 



To obtain a/ we use equation 1.78 from Rybicki & Lightman (1985) i.e., 

^® Notice that we are using just thermal broadening in the force multiplier calculations. We are not using turbulent velocity in these 
calculations. The use of turbulent velocities in this work is to simulate the effect of the acceleration between layers in the obtention 
of absorption profiles using our multilayered nearly relativistic outflow model ^ Moreover, by using only thermal broadening in our 
calculations, we can directly compare our results to those of other authors. 



22 




Fig. a. 4. — A comparison between our force multiplier calculations (solid lines) and those of Everett 2005 (dotted lines) and Arav et 
al. 1994 (dashed lines). The upper-left panel shows the continuum force multiplier as a function of the ionization parameter U . The other 
three panels show the line force multiplier as function of U for three different values of t. 
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Comparison with other results. 

The force multiplier mainly depends on the SED of the source, the composition of the gas, and its ionization 
parameter. In general a soft spectrum irradiating a gas that is nearly neutral (i.e. low ionization parameter) will 
strongly accelerate the gas and the force multiplier will be relatively large. 

When the gas is highly ionized it absorbs less radiation and is subject to a weaker radiative driving force (i.e., the 
force multiplier is relatively small). 

Figure A. 4 shows the estimated force multiplier assuming a Mathews-Ferland SED. Figure A. 4 has four panels, 
the upper-left panel represents the continuum component of the force multiplier, mainly produced by bound-bound 
transitions. The other three panels in Figure |A.4| show the contribution to the force multiplier from bound-free 
transitions for three different values of the dimensionless factor t; log t — —6, —7 and —8. The line force multipliers 
increase as t decreases. This increase is noticeable when the me dium is not highly i onized (log U < —1; log ^ < 0.3). 
Our calculations are in good agreement with those performed by Arav et al. (19941; Everett (2005). 



FORCE MULTIPLIERS CALCULATIONS OF MATERIAL ILLUMINATED BY ABSORBED POWER-LAWS AND 

MATHEWS-FERLAND SED 

A warm absorber shielding the outflow may increase the effectiveness of radiation driving. In this appendix we 
use the same unabsorbed SEDs used in ^ to estimate the amount of shielding required to significantly increase the 
radiative driving force. Additionally, we estimate the differences between the radiative forces produced by a soft 
incident spectrum and a hard incident spectrum for different ionization levels of the absorbing outflow. As in Q we 
assume that the warm absorber shielding the outflow has log (iVn/cm^^) = 23. 

In Figure A.5| we have plotted the bound- free, Mc, and the bound-bound, Ml, components of the force multiplier 
as a function of the ionization parameter for soft (F = 2.1) and hard (F = 1.7) absorbed power-law SEDs. In the 
upper-left, upper-right, lower-left and lower-right panel we have assu med absorbing shields with log {Nn/cm^^) = 23 
and log = 0.6, 1.0, 1.3, 1.6, 1.8 and 2.1, respectively. In Figure 1X3] we show that Ml is significantly larger for 
incident SEDs with soft spectra (i.e., F — 2.1) than for incident SEDs with hard spectra (i.e., F = 1.7) especially at 
low ionization levels of the shielding gas. For val ues of log $„v, 5-, 0-6 we do not expect strong radiative driving given 
the heavy attenuation of the continuum (see e.g., Chelouche fc Netzer||2003" ). We also note t hat t he force multipliers 
Mc and Ml are significantly reduced for low ionization levels of the shield as shown in Figure A.5[ Additionally, when 
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Fig. a. 5. — The bound-free, Mq, and bound-bound, Mi^, components of the force multiplier are shown as a function of the ionization 
parameter. Force multipliers are calculated for a power-law SEDs with photon indices of F = 2.1 (solid lines) and F = 1.7 (dotted lines). 
In the upper- left, upper-right, lower-left and lower-right panel we have assumed absorbing shield where the soft and hard SEDs have been 
attenuated by a warm-absorber shield with log (A^jj/cm^^) = 23 and log ^^.j^ = 0.6, 1.0, 1.3, 1.6, 1.8 and 2.1 respectively. 
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Fig. a. 6. — The bound-free, Mq, and bound-bound. Ml, components of the force multiplier are shown as a function of the ionization 
parameter. Force multipliers are calculated for Mathews-Ferland SEDs with X-ray photon indices of F = 2.1 (qqx = —1.52; solid lines) and 
F = 1.7 (qox = —1.41 dotted lines). In the upper-left, upper-right, lower-left and lower-right panels we have assumed that the SEDs have 
been attenuated by shields with log (TVjj/cm"^) = 23 and log ^^j^ = 1.0, 1.2, 1.9, 2.7, 3.1 and 3.8, respectively. The default parameters of 
a Mathews-Ferland SED are F = 1.7 and oiox = -1.41. 

log ^sIj > 2.0 the shielding becomes transparent to the radiation and therefore the force multiphers are similar to the 



ones in the unabsorbed case. Finally from Figure A. 5 we see that when the ionization parameter of the shield lies in 
the range of 1.0 < log ^^j^ ^1.8 the line force multiplier of the outflowing material is greater with the presence of an 
absorbing shield than without one. 

In Figure A. 6 we have plotted Mc, and Ml as a function of the ionization parameter for soft (F = 2.1; aox = —1.51) 
and hard (1' = 1.7; aox = —1.41) absorbed Mathews-Ferland SEDs. In the upper-left, upper-right, lower-left and 
lower-right panels we assume shields that attenuate the SEDs with log (A^n/cm^^) = 23 and log ^^jj = 1.0, 1.2, 1.9, 
2.7, 3.1 and 3.8, respectively. In the case of a Mathews-Ferland SED incident on a shield with an ionization parameter 
in the range of 1.2 < log ^^j^ < 3.8 we find that the line force multiplier Ml for soft SEDs (i.e., F = 2.1) extends to 
larger value s of the ionization parameter of the outflowing material compared to hard SEDs (i.e., F = 1.7). 

In Figure [A77 we have plotted Mc, and Ml as a function of the ionization parameter of the outflow for soft (F — 1.7; 
ttox = —1.51) and hard (F = 1.7; Ofox = —1.41) absorbed modified Mathews-Ferland SEDs. In the upper-left, upper- 
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Fig. a. 7. — The bound-free, Mq, and bound-bound, Mi^, components of the force multiplier are shown as a function of the ionization 
parameter. Force multipliers are calculated for Mathews-Ferland SEDs with X-ray photon indices of ctox = —1.52 (F = 1.7; solid lines) 
and Oox = —1.41 (F = 1.7; dotted lines). In the upper-left, upper-right, lower-left and lower-right panels we have assumed that the SEDs 
have been attenuated by absorbing shields with log (A'^jj/cm"^) = 23 and log ^^jj = 1.0, 1.2, 1.9, 2.7, 3.1 and 3.8, respectively. The default 
parameters of a Mathews-Ferland SED are F = 1.7 and Oox = —1.41. 
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Fig. a. 8. — Temperature of a warm absorber (shield) as a function of its ionization parameter (illuminated side) for a warm absorber 
with log {NYi/cm~^) = 23. In the left panel the illuminating source is a power law with F = 2.1 (solid line) and F = 1.7 (dotted line). 
In the right panel the illuminating source is a modified Mathews-Ferland SED with F = 2.1 and aox = —1.52 (solid line), F = 1.7 and 
Qox = —1.41 (dotted line), and F = 1.7 and Oox = —1.52 (dashed line). 

right, lower-left and lower-right panels we assume shields that attenuate the SEDs with log (A^n/cni^^) = 23 and 
log = 1.0, 1.2, 1.9, 2.7, 3.1 and 3.8, respectively. In the case of a modified Mathews-Ferland SED incident on a 
shield with an ionization parameter in the range of 1.0 < log < 3.9, we find that the line force multiplier is 
significantly larger for soft SEDs (i.e., ckox — —1.52) than for hard SEDs (i.e., — —1.41). We note that the line force 
multiplier for the case of a soft Mathews-Ferland SED incident on a shield with an ionization parameter of log ^gij=1.9 
has a value larger than ^1000 for ionization parameters of the outflowing absorber in the range of < log ^ < 5. For 
a shield with log ^gj^ < 1 we predict that the continuum becomes too attenuated to allow effective radiative driving. 
When log ^^jj > 3.1 the shield becomes transparent to the incident radiation and ther efore the forc e multipliers are 
similar to th e on es estimated for the unabsorbed case (see lower right panels of Figures A. 6 and A.7). 

In Figure A. 8 we plot the temperature of an absorbing shield as a function of its ionization parameter (on its 
illuminated side) for a shield with log (iVn/cm^^) — 23. From the analysis in the last three paragraphs and from 



Figure A. 8 we conclude that the most powerful winds arc probably produced at temperatures in the range 4.5 < 



log T < 5.2. This occurs for power-law SEDs incident on shields with 1.0 < log ^gjj < 1.8 and for Mathews-Ferland 
SEDs incident on shields with 1.2 < log ^ 3.1. 



These are temperatures at the illuminated face of the shield (output from cloudy). 



